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Introduction. 

§ 1. The object of the present paper can be best explained by referring to a sentence 
which occurs in a paper by Professor G. H. Darwin. ^^ This is as follows :-— 

^' The principal question involved in the nebular hypothesis seems to be the 
stability of a rotating mass of gas ; but, unfortunately, this has remained up to now 
an untouched field of mathematical research. We can only judge of probable results 
from the investigations which have been made concerning the stability of a rotating 
mass of liquid.'' 

In so far as the two cases are parallel, the argument by analogy will, of course, be 
valid enough, but the compressibility of a gas makes possible in the gaseous nebula a 
Avhole series of vibrations which have no counterpart in a liquid, and no inference as 
to the stability of these motions can be drawn from an examination of the behaviour 
of a liquid. Thus, although there will be unstable vibrations in a rotating mass of 
gas similar to those which are known to exist in a rotating liquid, it does not at all 
follow that a rotating gas will become unstable, in the first place, through vibrations 
which have a counterpart in a rotating liquid : it is at any rate conceivable that the 
vibrations through which the gas first becomes unstable are vibrations in which the 
compressibility of the gas plays so prominent a part, that no vibration of the kind 
can occur in a liquid. If this is so, the conditions of the fox^mation of planetary 
systems will be widely differeni: in the two cases. 

With a view to answering the questions suggested by this argument, the present 
paper attempts to examine in a direct manner the stability of a mass of gravitating 
gas, and it will be found that, on the whole, the results are not such as could have 
been predicted by analogy from the results in the case of a gravitating liquid. The 

"^ "On the Meclianical Conditions of a Swarm of Meteorites, and on Theories of Cosmogony," 'Phil. 
Trans.,' A, vol. 180, p. 1 (1888). 
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2 MR. J. H. JEANS ON THE STABILITY OF A SPHEEICAL NEBULA, 

main point of difference between the two cases can be seen^ almost without 
mathematical analysis, as follows : — - 

§ 2. Speaking somewhat loosely, the stability or instability may be measured by 
the resultant of several factors. In the case of an incompressible liquid we may say 
that gravitation tends to stability, and rotation to instability ; the liquid becomes 
unstable as soon as the second factor preponderates over the first. The gravitational 
tendency to stability arises in this case from the surface inequalities caused by the 
displacement : matter is moved from a place of higher potential to a place of lower 
potential, and in this way the gravitational potential energy is increased. As soon 
as we pass to the consideration of a compressible gas the case is entirely different. 

Suppose, to take the simplest case, that we are dealing with a single shell of 
gravitating gas, bounded by spheres of radii r and r + dr, and initially in equilibrium 
under its own gravitation, at a uniform density pQ. 

Suppose, now, that this gas is caused to undergo a tangential compression or 
dilatation, such that the density is changed from 

Pi) to pq + ^p&i> 

where pu is a small quantity, and S,, is a spherical surface harmonic of order tl 

It wdll readily be verified that there is a decrease in the gravitational energy of 

47rr^ {chf t .^'^ f [ S/ sin 6 dd d<l>. 

\ ^ ' t* "i J. / J J 

As this is essentially a positive quantity, we see that any tangential displacement 
of a single shell will decrease the gravitational energy. 

This example is sufficient to show that when the gas is compressible, the tendency 
of gravitation may be towards instability. The gravitation of the surface inequalities 
will as before tend towards stabiHty, but when we are dealing with a gaseous nebula, 
it is impossible to suppose that a discontinuity of density can occur such as would be 
necessary if this tendency were to come into operation. Rotation as before will tend 
to instability^ and the factor which makes for stability will be the elasticity of 
uUe gas. 

We can now see that there is nothing inherently impossible, or even improbable, 
in the supposition that for a gaseous nebula the symmetrical configuration may 
become unstable even in the absence of rotation. The question which we shall 
primarily attempt to answer is, whether or not this is, in point of fact, a possible 
occurrence, and if so, under what circumstances it will take place. To investigate 
this problem, it will be sufficient to consider the vibrations of a non-rotating nebula 
about a configuration of spherical symmetry, 

I 3. Unfortunately, the stability of a gaseous nebula of finite size is not a subject 
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which lends itself well to mathematical treatment. The principal difficulty lies in 
finding a system which shall satisfy the ordinarily assumed gas equations, and shall 
at the same time give an adequate representation of the primitive nebula of 
astronomy. 

If we begin by supposing a nebula to consist of a gas which satisfies at every 
point the ordinarily assumed gas equations, and to be free from the influence of all 
external forces, then the only configuration of equilibrium is one which extends to an 
infinite distance, and is such that the nebula contains an infinite mass of gas. The 
only alternative is to suppose the gas to be totally devoid of thermal conductivity, 
and in this case there is an equilibrium configuration which is of finite size and 
involves only a finite mass of gas. But the assumption that a gas may be treated as 
non-conducting finds no justification in nature. When we are dealing, as in the 
present case, with changes extending through the course of thousands of years, we 
cannot suppose the gas to be such a bad conductor of heat, that any configuration, 
other than one of thermal equilibrium, may be regarded as permanent. 

Professor Darwin has pointed out that a nebula which consists of a swarm of 
meteorites may, under certain limitations, be treated as a gas of which the meteorites 
are the ^' molecules."^ In this quasi-gas the mean time of describing a free path must 
be measured in days, rather than (as in the case of an actual gas) in units of 
10""^ second. The process of equalisation of temperature will therefore be much 
slower than in the case of an actual gas, and it is possible that the conduction of heat 
may be so slow that it would be legitimate to regard adiabatic equilibrium as 
permanent, t 

Except for this the mathematical conditions are identical, whether we assume the 
gaseous or meteoritic hypothesis. The present paper deals primarily with a nebula in 
which the equilibrium is conductive, but it will be found possible from the results 
arrived at, to obtain some insight into the behaviour of a nebula in which the 
equilibrium is partially or wholly convective. 

S 4. Whether we suppose the thermal equilibrium of the 2:as to be conductive or 
adiabatic, we are still met by the difficulty that the gas equations break down over 
the outermost part of the nebula, through the density not being sufficiently great to 
warrant the statistical methods of the kinetic theory. This difficulty could be avoided 
by supposing that the nebula is of finite size, and that equilibrium is maintained by 
a constant pressure applied to the outer surface of the nebuhi. If this pressure is so 
great that the density of gas at the outer surface of the nebula is sufficiently large to 
justify us in supposing that the gas equations are satisfied everywhere inside this 
surface, then the difficulty in question will have been removed. On the other hand, 
this pressure can only be produced in nature by the impact of matter, this matter 

■^ G. H. Darwin, he, eii,^ ante. 
t Ibid,, p. 64. 
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C3nslstliii>' either of molecules or meteorites, so that we are now called upon to take 
account of the gravitational forces exerted upon the nebula by this matter. This 
whole question is, however, deferred until a later stage ; for the present we turn to 
the purely mathematical problem of finding the vibrations of a mass of gas which is in 
equilibrium in a spherical configuration. We shall consider two distinct cases. In the 
first, equilibrium is maintained by a constant pressure applied to the outer surface of 
the nebula, this surface being of radius K^. In the second, the nebula extends to infinity, 
and it is assumed that the ordinary gas equations are satisfied without limitation. We 
suppose for the present that the gas is in thermal equilibrium throughout. It is not, 
however, supposed that the gas is all at the same temperature ; to make the question 
more general, and to give a closer resemblance to the state of things which may be 
supposed to exist in nature, it will be supposed that the gas is collected round a solid 
spherical core of radius Rq, and the temperature will be supposed to fall off as we 
recede from this core to the surface, the equation of conduction of heat being satisfied 
at every point. We shall also suppose that the gas is acted upon by an external 
system of forces, this system being, like the nebula, spherically symmetrical. The 
reason for these generalisations will be seen later ; it will at any time be possible to 
pass to less general cases. 



The Criterion of Stability. 

The Principal Vibrations of a Spherical Nehida. 

§ 5. We shall take the point about which the nebula is symmetrical as origin. It 
will be convenient to use rectangular co-ordinates x, y, 7i^ in conjunction with polar 
co-ordinates r, ^, (^. We shall imagine the nebula to undergo a small continuous 
displacement ; let the components of this be ^, tj, ^, when referred to rectangular 
co-ordinates, and u^ rv^ riv sin when referred to polars. Thus the point initially at 



x^ 2/5 z or r^ 0^ (j) 



is found after displacement at 



^ + ^, y + V, ^ + C or r + 11, + V, <j> + u\ 
The cubical dilatation of this displacement will be denoted by A, so that 

dx dy dz 



7^ |: (^'«) + ,i^ fe (^' «i^ ^) + f^ 
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In general we shall denote the density by />, pressure by ot, temperature by T, 
total potential by V, coefficient of coyid action of heat by ac, and the gas constant by X, 
the last of these being given by the equation 

t3- = XTp (1). 

In the equilibrium configuration each of the quantities just defined is a function of 
r only. 

If c is any one of these quantities, we shall denote the 

Value of c in the equilibrium configuration, evaluated at x, y^ z^ by Cq. 

Q-isijiaceoL jj ,j jj jj Cr\ ~p" c . 

at x + ^, y + 77, ;^ + ^ by Co + Cp 



?? 55 



;? J? J5 5? ?? 



The quantities Gq, c , c^ are, of course, not independent. Since Cq + c^ is the same 
function of cc + ^5 ^ + '*?? ^-^ + C? ^s is Cq + c' of x, y^ 7^ we have, as far as the first 
order of small quantities, 

Co + e^ = .0 + c + ^9^ + ^ 9^ + rf , 



or, since Cq is a function of r only, 



U J^ — — I./ j (.</ , • « « « » . , , , , \Zj J , 



§ 6. From the equation of continuity we have at once 

Pi — *" """ Po ^ ' * • • • * * • • • • V*-^/* 

Since X remains the same throughout the motion of any given element of the gas, 
Hence, from equation (l). 



^0 + ^1 = ^0 (To + Ti) (po + Pi): 



giving as the value of txT^ 



(37 



1 == Xq (TjPo + TqP^) ^ K,po(^i — ATq) ...... (5). 



So long as we confine our attention to a single, element of the gas, the coefficient of 
conduction of heat is proportional to the square root of the temperature, and is 
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independent of the density. "^^ We therefore have, as far as the first order of small 
quantities, 






Lastly Y\ regarded as the difference between Vq + V^ and Vq, is seen to be the 
potential of a volume-distribution of matter of density p\ to which must be added : 

(i.) The potential of a surface-distribution over the sphere r = E-q, the surface 
density being 



Mpo — ^o)." 



r — Eo' 



where o-q is the mean density of the core, and 

(ii.) The potential of a surface-distribution over the sphere r 
density being 



Rt, the surface 



where cr^ is the density of the medium (if any) outside the nebula. 

§ 7. We are now in a position to handle the equations of motion, and of conduction 
of heat. For the element which, in the undisturbed state, is at x^ y^ z^ the equations 
of motion are three of the type 



3^? 9 

8^3 3^ ^ ^ ^ (Po + pO 



. g^K+t^O 



(7) 

» » e » \ w f * 



Transforming to polar co-ordinates, these equations are equivalent to 



r 



3^'V 
3? 



1 3V' 

r do 



1 dm^ 
p^T do 



B . 



(8). 



. . V /• 



. „dhv 1 3V' 1 Bey' 

dt^ T sin d<p p^r sin 3^ 



As an equation of equilibrium, we have 



dYr 



1 3 



m 







Po 9^ 







and with the help of this, equation (8) reduces to 



3^ 



•It 



dp 



3Y' 1 dm^ p' dm^ 

dr Pq dr p^" dr ' 



as far as the first order of small quantities. 



. . (10). 



8 O O « C 



(11): 



« d 9 b « 9 IJL dLJ /* 



^ BoLTZMANN, * Vorlesuugen iiber Gastheorie/ vol. 1, § 13. 
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Let us write 



so that 



X 



v — 



Z*5" 



Po 



• • « 



• • • • 



9% 


dY' 


1 vm^ m^ upQ 

Po 3'^^ Pi)^ 3^' ' 


or 


dr 



then equation (12) becomes 



"■ 2 —' X^ 'r- 'oAP ;^,; — ^' 



dt^ dr 



Po 



dr ' ' 



(13): 



and, by the use of equation (2), this is seen to be equivalent to 



3% ^X i_ 1 / 

dt^ dr Pq y 



dm 







— 'CiTi 



dr I ' 



' (14) 



Equations (9) and (10) now take the simple forms. 



d^v 
dt^ 






33, 



to 



1 



^X 



dt^ r^ sin^ 6 3^ 



From these last two equations, we obtain at once 



:^2 



dt^ iskr^al(''^'''^^ + 3^. 



% 



1 3 / . ^ a^\ 1 33 

r^ sin ^ dd \ dO j r^ sin^ dcj)^ 



or, what is the same thing, 



d^ ( 1:^ 



: (^^^^0 \ = "^^x 



1 A /..2 3%' 

r^ dr 



• • » 



. (15). 



§ 8. The equation of conduction of heat is, as far as the first order of small 
quantities, 



dt 



JltJL r\ . I ^^P ' 






1 
p 



J 3 / 8T\ , a / aT\ a / aT\| 

_3^\ dxj dy\ dijl ^ dz\ dz I \' ' ' ^ ^' 



in which p, k, T stand for Pq + p\ ^o + ^^ '^o + '^^ respectively. The notation is that 
of KiRCHHOFF ; the equation may either be written down from first principles, or 
regarded as a simplified form of Kirchhoff's general equation/^' 
Since there is thermal equilibrium in the undisturbed configuration, 



a/ 9To\ a / siv , 



dx \ '^^ dx 



9 ' 3n\ _ 



d% \ '^^ d% 



» « <! 



• • (17). 



* KiRCHHOFF, ' Vorlesungen fiber die Theorie der Warme, p. 118. 
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Hence equation (16) reduces to the form 






1 r a / aT'\ , a 

— — J — . j 1^ , , j I — , 

Po ISb\ ° a»/ a?/ 



aT'^ 



a 



!^T/ 



"0 ay i + a 






K 



"-p 7 



8 



K 



; 3X 







dw \ 3 



\ 



te 






d 



/C 



9^ 

^ Ox A 



4- 



a 



K'?^^U 



dy \ dy j dz \ c) 



vz 



. (18). 



Since Kg, Tq are functions of r only, the bracket on the right-hand side of this last 
equation again reduces to 



a«n aT' 



&!• or ^ dr 



a«' aT 



or 



4" K V^T 







9 s » 



s # • « \/5 



and, cleared of accented symbols by the use of equation (2), this takes the form 



d?' or ^' ^ err or ^ 







3r \ 8?" 8r 



8'?^ J 8/^:^) 8T0 
or [ or or 






/8T. 



- - 



8/c 



\^ 8/" 










Orn 



4_. I/. 







if I"" 

01 



t-70 



-^0 'a, v~ t. 



* « 



Now equation (17) can be written in the form 



(20). 



okq BTq 



dr dr ^ ^^ V ±0 »^ U 



« e s 9 e 5 d 



(21): 



whence, by differentiation with respect to r. 



or V <^/^ o''" I or ^ or ^^ 



• a 



« » • \ <^ JmJ f 9 



With the help of equation (22), the bracket in the second hne of (20) reduces to 






"'■^ or 



r 



fr 



while, with the help of (21), that in the third line becomes 



/j Kr\ O JL r\ 



r c 



dr 



Again, if we substitute for /Cj the value found for it in equation (6), the two last 
terms in the first line of (20) can be transformed as follows : 



^ ^_o J_ ,^ X7-2T ,^ '^ / ^l \ ^-^^o _i tl 

),, 97 + '^iV lo - '^0 37, i 2T„; "dr + 2T„ 



01 



0/ 






a 



r ar +''oV-Jo 



>■ ' 



and the last bracket vanishes by equation (21). 
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Collecting results, and substituting for p^ from ec^uation (3), we find that equation 
(18) takes the form 



■^^Po ~i{f ~r ^i" a7 



1 
P 



B/^o aTj 



a /Ti\aT( 



or 



4 dib 



3To [2% 



IT \ 9 






. (23).- 



§ 9. In addition to the volume-equations whicli have just been found, there are 
certain boundary conditions which must be satisfied. These are as follows : 

(i.) The pressure must remain constant at the outer surface, so that we must 
liave 

(ii.) The temperature must remain unaltered at r = Eq, or else the flow of 
temperature across the surface 7' — Rq must remain 7iiL These two suppositions 
require respectively 



[Tl]r=Ro = 0^ ^1^ 



dr 



0. 



u'-Ro 



(iii.) A similar temperature condition must be satisfied at 7* =: Rp 

(iv.) The kinematical and dynamical boundary conditions at the sui'face 7' = Rq 
must be satisfied. These express that the normal velocities shall be continuous at 
this surface, and that the motion of the rigid core shall be such as would be caused 
by the forces acting upon it from the gas. 

§ 10. Equations (14), (15) and (23) give the rates of change in u, A and T\ in terms 
of these quantities. Hence these equations enable us theoretically to trace the 
changes in u, A and T^, starting from any arbitrary values of u^ A, T^, du/dt and 
dA/dt, which are such as to satisfy the boundary conditions. 

Imagine initial values of ti, A, T^, du/dt and dA/dt, in which the latitude and 
longitude enter only through the factor S;„ where S,^ is any spherical harmonic of 
order 7i. Then it can be shown that the solution through all time (so long as the 
squares of the displacement may be neglected) is such that the latitude and longitude 
enter only through the factor S,,. For, assuming a solution of this form, the value of 
V found in § 6 will contain S;, as a factor, as will also p^, -^i, m (equations 3, 5, 2) 
and X (equation 13). The same is true of V% V^Tj^ and V^i/, since 



V^ [/(r) SJ 



^ ^ f..,df(r)\ n(n+l)f('r) 



r^ dr 



dr 



o 



s 



Uf 



* Q- 



Sections 5-8 were re- written in November, 1901. I take this oj)portunity of expressing nw thanks 
to the referee for the care and trouble which he h^is bestowed upon my paper. To him I am indebted for 
several improvements in these four sections, in particular for the present form of equation. (23), and also 
for the removal of a serious inaccuracy from my original equations. 

VOL. OXCIX.-^^ -^A. 
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where y(r) is any function of r. It therefore appears that every term in equations (14), 
(15) and (23) will contain S^^ as a factor. Dividing out by this factor, we are left with 
equa^tions which do not involve 9 and ^, and this verifies our statement. 

§ 11. It therefore follows that there are principal vibrations'^^ in which ii, A and T^ 
are of the form 



A = BS,;,e^^^ 



T, =^ CS,e^ 



. (25), 

• (26), 



lu which A, B, are functions of r only. The relations between A, B, C and p must 
be found from the equations (14), (15), (23), and the boundary conditions. 
The value of p for the vibration just specified is 



Pi 



li 



^0 

dr 



Apo + u 



dpQ 







dp 



A ^ + Bpo j S>^. 



We shall in future drop all zero suffixes^ there being no longer any danger of 
confusion. Calculating V after the manner explained in § 6, we find {of, Thomson 
and Tait, 'Nat. Phil./ § 542), 

where 



V 



47r 



{in + l)r^^+i 



Rn 



A ^ + Bp ) r«+^ c^r 



:A(p--o-o)r"+^^J,.i,„ 



+ /7 



(271 + 1) 



{-.C(^:f+»'');^-+ 






r=lli J 



. . . (27). 



We liave further^ by equations (2) and (5), 



d 



m 



ZtTi 



It 



m 



dm 



dr 



Xp (C ~ BT) S„e'>' - A -^ S«e'^^ 



and hence we obtain (equation 13) 



Vvliei'e 



X = 'FS„0^", 



F = V - X (0 



rt^ 



) 






A du 



d 



m 



I' 



i > 



(28). 



Substituting the assumed solutions for u, A and T^, and the corresponding values 
for x? Pi^ ^1^ 1^^ equations (14), (15) and (23), and dividing throughout by the factor 
B^f,e'P\ we find the relations 

* In order to avoid circumlocution, we shall find it convenient to use the terms '' principal co-ordinate " 
and '^ principal vibration," although we are ignorant as to whether the nebula is stable or unstable. It 
will ultimately be found that we only apply our results to nebulae which are either stable or in the limiting 
state of neutral equilibrium* 
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p~A 



d¥ 
dr 



p^B 



B dm 

p dr 

7"^ ar ^ ^ 



X dp 
p dr 



(C-BT) ..... (29)^ 



n(n + 1) 






F 



^>p{MpB + C„C) = 



d/c dG 
dr dr 



+ 



"{f('^f)-"(»+')° 



K 



dT 

dr 



1 d 



r^ dr 



/ 

o dK\ 
dr J 



wm^imm* §/" 



d f C \ ^T 



4 dA. 
r dr 



dr \ 2T/ dr 



The boundary-equations found in § 9 reduce to the following : 



(i) [c - bt; 



■=E, 



. . . 



• « • 



• « « • « 



• i 



(ii) C,^R„ = or 



'dXS 
d>T 






r=Ko 



(iii) Equations similai" to (33) at r =: R, . . . . , 



... 



(30), 



(31). 



. . . (32), 



. (33), 



• . ♦ • \^O^j^ 



(iv) {A)^^R^ = 0, when n is different from unity^ or a more complex equation in the 

VCLk50 \j\. iv •-»»-• J. t , , .-» . • .. a a * « . . . « » « \ tj O ) , 



§ 12. From the manner in which the analysis has been conducted^ it will be clear 
that every principal vibration must either be one of the class just investigated, or 
else a vibration such that it, A, and T vanish everywhere. 

For the latter class of vibration there are no forces of restitution. Thus the 
frequency of vibration is zero, and the motion consists of the flow of the gas in closed 
circuits, this flow being entirely tangential, and the gas behaving like an incom- 
pressible fluid. Obviously these steady currents are of no importance in connection 
with the question of stability or instability. 



Discussion of the Frequency Equation. 

§ 13. Returning to the class of vibrations in which u, A, and T do not all vanish, 
we have seen that the frequency equation is found by the elimination of F, A, B, and 
C from equations (28) to (35). Now p only enters into three of these equations : 
namely (31), in which it enters through the factor ^p, and (29) and (30), in which it 
enters through the factor -— p^ or {ip)^. Regarding ip^ A, B, C, and F as unknowns, 
it will be seen that the coefiicients which occur in equations (28) to (35) are all real. 
The four volume equations enable us to determine A, B, C, and F except for certain 

c 2 
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constants of integration^ and tlie values of these quantities will be wholly real if ip is 
real. The boundary-equations enable us to determine the constants of integration 
and also provide an equation for ip. Every term in these equations will be real if ip 
is real. Hence the frequency equation can be written in the form. 

fiip) = 0, 

where f{x) is a function of x in which all the coefficients are real, these coefficients 
being functions of n and of the quantities which determine the equilibrium configura- 
tion of the nebula. 

It follows that the complex roots of ^p will occur in pairs of the form 

ip — 7 ± ^' 8. 

where y and S are both real. There ma)?- also be roots for which ip is purely real^ so 
that S = Oj and y exists alone. 

The vibration corresponding to any root is stable or unstable according as y is 
negative or positive. 

If the equilibrium configuration of the nebula changes in any continuous manner, 
so as always to remain an equilibrium configuratioUj the values of ip will also change 
in a continuous manner, and for physical reasons these values can never become 
infinite. Hence, if the configuration of the nebula changes from one of stability to 
one of instability, it must do so by passing through a configuration in which there is 
a vibration for which y =: 0. 

§ 14. For the present we shall not discuss the actual stability or instability of any 
configuration, but shall examine under w^hat circumstances a transition from stability 
to instability can occur. 

We therefore proceed to search for configurations in which there are vibrations 
such that y =: 0. Now for such a vibration we have either a root of the frequency 
equation ^9 = 0, or else a pair of roots of the form ip ^^ ^ 2'S. 

In the latter case the corresponding vibration is one in which a dissipation of energy 
does not occur. A necessary condition for such a vibration is that no conduction of 
heat shall take place. Hence both sides of the equation of conduction of heat 
(equation 31) must vanish. Excluding adiabatic motion (represented by the 
vanishing of the factor MpB -f C„C), this condition, compels us to take 



fj •~-»« \j 



together with 



d. Jr + 7^ i*- r ^'^^ - « (^ + 1) C j + K ^^, {^^^ -^_ 

fZT r 1 d I ^ dk\ 4 dK n {n + 1) - 2 , ' 

dr 'fi clr V dr / r dr r^ 



. . (36). 
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la 



Thus vibrations for which y = 0, if they exist, must satisfy equations (32) to (36), 
and also equations (29) and (30), in which jd is put equal to zero, and equation (28), 

The case of n = will be considered later (§ 28). Excluding this for the present, 
we find that puttings = in (30) leads to 



F=:0 



♦ «»««» 



Equation (29) now reduces to 



« « » t « \tjfif 



B -^'^ 4- X ^? (C ~ BT) 

rim ri V ^ ' 



dr 



d}'- 



, 



«««««« 



. (38), 



or, replacing ot by its value XTp, 



Bp 'I {XT) + X '^'J C := . 



dr 



dr 



« * » « f 3 t 



(89), 



Equation (28) becomes 



V = X (C - BT) 



k cl 



m 



dr 







* « » « 



« « 



(40), 



atid tbe (iliiTiination of Q — BT from tbis equation and (38) leads to tbe equation 



dp^. 



dr 



Y 



\ 



dp 

dr 



A ~ -f- Bp 



1 dm 
p dr 



e e « • 9 



Substituting for V from equation (27), this becomes 



r 



4:7r 
(2n + l)r*H-i ^j^ 



M 



4.7rf 



11 



(271 + 1) 



Aj|- + B/3)r'"+3c/r4 



■'' A ^ + Bp^ ~ '- 

dr ^ / r'"^ 



A (p — o-q) r'''^' 



'A (p - a J 



2 



1 



r=B. J 



,wa—\ 



'==R, 



A '^- + Bp V^ -~ Z^'' 

r/?' '^ } p dr I dr 



• ■ • « 



» a « » 



(41). 



(42). 



§ 15. With a view to transforming this equation, let us consider the equation. 



47r 



{2n + l)7-+i LJb, 



!'■ Jr»+2,?r + Kol + 



^nrf 



It 



(2n + 1) 



•Ri T 



Ij (4o). 



in which J and L are any functions of i% and Kq, K^ are constants. If we multiply 
by v'''^^, and differentiate with respect to r, we obtain, after some simplification, 



4777 



,%t 



<^ 



Ri J 



d 



.. '? 



= X(^^"'') • . . . • . (44), 



while by multiplying (43) by r'^'' and differentiating, we obtain in a similar way 



4t.lT 



"T 



^ I JR, 



Jr«+2 dr + K 



fZ 







Q 



«!? 



: (Lr-«) (45). 
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Divide (44) by r^^' and differentiate with respect to r, then 






or, writing ^ for Lr, and simplifying 



#1 fh (n 4- 1) 



dr^ 






i 



4#rJ ♦ . (47)j 



and this same equation conld hare been deduced from. (45) instead of (44). 

Equation (47) is more general than (43) since the two constants K^^ Kj have 
disappeared. In fact equation (47), being a differential equation of the second order, 
will contain two arbitrary constants in its solution, and these correspond to the two 

missin^t constants Ka and Ki, We can, however, determine Ka, K-, in terms of 



these two arbitrary constants, and if these constants are chosen so as to give the 
right values for Kq, Kj, the solution of (47) will be equivalent to the original 
equation (43). 

To determine Kq^ K^, put r =: Ej in (44) and we obtain 



4:frK,' 



1 d 



« * » * 



'*^=B 



♦ (48)j 



and similarly from (45) 



h 



4YT 







7 



Mn~h% 



Al 



zii''~ ■'''') 



di 






(49). 



Hence we see that equation (43) is exactly equivalent to the three equations 
(47)5 (48), and (49), 

I 16. Comparing (42) with (48), it appears that (42) is exactly equivalent to the 
following equations :- — 

(50). 



> i i dp , ^ \ r dm I dp 

1^ = A -r + Bp - ~ / -f- 

^^''" ' I p dr / m' 



di 



«#:»«»£ 



#f n (b +1) 
dr^ 



n-in 



i 



/ 7 \ 

iwr A 'v + Bp 

\ dr ' 



/ 



s ♦ i «J J. L 



4^J\.| 



X (.V / g» \ 






^ciii f~i^ 



■f=B, 



"A(p-o-J 



/jjji 1 



r=Ri 



iijZija 



4irK. 



•0 



d 



1 1' / 



r=Rf 



The right-hand member of (51) is equal to 



A (p -«- (Tq) f 



, J, dp 11 dm 
dr I p dr 



tfi + s 



f=:R, 



I Do j» 
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so that if we introduce a new quantity u, defined by 

dp I dm > 

ic^27rpr^ j^Iy: (54), 



dr / dr 



equation (51) may be written in the form 



r^^= {n{n + 1) — 2u} ^. ,.,.,, , (55). 



The sohition of tliis will be of the form 

I = Ei^i (r) + Eoc/>^ (f ) ........ (56), 

in which Ej, E^ are constants of integration. We have, from the definition of |, 

A'^^ + Bp=z^^^^^= . , . . (57), 

and the elimination of B from this equation and (39) gives 

^^"C = jr(XT)|A^5-£,,{EAW + E,<^30-)}| • . . (58). 



dr dr 



dr 2wi' 



If we imagine this value for C substituted in equation (36), we shall have a 
difterential equation of the second order for A. The solution of this will be of the 
form 

A-EJ,(r) + E,/,(r) + Ej;(r) + E,./l(r) .... (59), 

in which E^ ctnd E^. are the new constants of integration. From this value of A we 
can doduce the values of B and C (equations (57) and (58)) without introducing any 
further constants of integration. 

Turning to the boundary conditions, we now find that tliere are six boundary- 
equations to be satisfied (equations (32), (33), (34), (35), (52), (53)) and only three 
arbitrary constants at our disposal, namely, the ratios of the four E s. If we 
eliminate these E's we shall be left with three equations to determine the configura- 
tion of the nebula at which instability sets in, and these equations will in general' be 
inconsistent. 

I 17. In order to put the right interpretation upon this result, it will be necessary 
to return to the general equations of free vibrations found in § 12. 

If we eliminate F from equations (29) and (30), we obtain 
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while equation (30) may^ with the help of (28), be written in the form 



V«-~. fi/f}^ /fill 



in which £ is now defined by 



=: X(C — BT) — ^- .. - + _Jl.^--^-_-^ J f^ «_ i.^^.3_^^ I , , . (62). 

/ p f^r ^i (^'^i -j- 1 J j^ (It J 

Substitnting for V from equation (27), and treating the equation so formed in the 
manner explained in § 15, w^e find, as the equivalent of equation (61), 

(i.) A volume equation^ analogous in form to (51), namely, 

#1 11 {n + 1) , / 



(I't 






4- == ^ 4^ri A ^-i^; + Bpj . , . . . (63), 



(ii.) Two boundary equations analogous in form to (52) and (53), . . (64)^ (65)» 



Thus the equations found in § 1 1 may be replaced by 

(a) Three volume equations, namely, equations (60), (63), and (31). 

0) Six boimdary equations, namely, equations (32), (33), (34), (35), (64), (65). 

We may conduct the elimination of B and C from the three equations (a) in a 
symbolic manner as follows : — 

Let D,, be a symbol ivhich is used to denote any Ihiear difierential operator of 
order n, the differentiations being wnth respect to n The symbol has reference 
solely to the order of the highest difierential coeflicient which occurs, and must in no 
case have reference to any particular differential operator. Thus we write D,, 
indiscriminately for every operator of the form 

,/« h') a,.. + fn-i h') g^:«::i + ^ ■. ^ 

The law.^s governing the manipulation of this symboFare as follow'^s : 

(ii) ■ D,,</j +" D,^c/> = !),//> (w/ > rn)'. 

(iii.) I);^ (I)^;^c/j) — \),f^^,^^i|y, 
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It must be particularly noticed that in general 

L/.;^(p ■— Jij-fid) ^^^^ L.)^(b, 

Corresponding, however, to any two specified operators of order n, say (D;,)|- and 
(D,,)2j it will always be possible to find two functions of r, say a and b, such that 

a (Dn)i <l> — h (Dn).2 4* = D./^_j«^ ....... (66). 

In terms of this operator, the three equations (a) (p. 16) may be written in the 

following forms : 

p^ (DgA + D]^B) + D()B + DqC = . (67), 

D,A + D,B + D,C + p' (D3A + D,B) - (68), 

ip yUgij -}- UqO) -f" DgO -f" DgA = . (69). 

Now D,, is commutative with regard to functions of r, and is of course commutative 
with regard to p. This enables us to eliminate B and C from the above equations. 
To make this clearer, consider a simple case, say the pair of equations 

DA _ T) 13 (^()\ 



D,A=p^D„B (71). 

If we operate on (71) with d/dr^ we get an equation of the form 

and from this and equation (70), we can, with the help ot the property expressed in 
equation {QQ), deduce an equation of the form 

From this and equation (71) we can in a similar way obtain an equation of the form 

T) A — Tl 1^ 4. ^2T) Jg 

We may regard this as an equation giving A, and substitute for A in (71). In this 
way we obtain 

and the elimination of A has been effected. 

It will be clear that throup^hout this elimination we have followed a method which 
would have been successful in eliminating A if d/dr had been regarded as a mere 
multiplier. The result of the elimination is accordingly exactly the same as might 
have been obtained directly from the original equations (70) and (71), by regarding 
the D's as multipliers and eliminating according to the ordinary laws of algebra, 

VOL. OXCIX. — A. D 
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It will now be apparent that we can eliminate any two of the three unknowns, A, 
B, and C, from equations (67)-{69) by this method. The difierential equation satisfied 
by the remaining unknown (say A) will be 



where, symbolically. 



D2, 



3? 



p^Di + Do, 
p^Dj + Dg, 



(73). 







D 



a 



ij^Dq ■+ D- 



74 



We may expand this determinant according to the rules ah*eady laid down for the 
manipulation of the D's, and so obtain 



A = '^j^^D^ + P^Dg + ^p^D^ + P^Dg + ipD^ + D 



4 



(75). 



§ 18, We can now see the explanation of the difficulty which occurred in § 1 
The occurrence of the term Dg in A points to a differential equation of the sixth 
order, which is satisfied by any one of the quantities A, B, or C in the general case, 
in which j) does not vanish. As soon, however, as p is put equal to zero, the 
expression for A reduces to D4^, and the differential equation is one of the fourth order 
only. It therefore appears that by putting p = before solving the differential 
equations, the order of these equations is reduced automatically, and two solutions 
are entirely lost from sight. 

These two last solutions, it is easy to see, are solutions which do not approximate 
to a definite limit, when p approximates to zero. The remaining four solutions will 
approximate to the same forms as would be obtained by putting p = before solving 
the differential equations. Thus, instead of equation (59), we must write the complete 
limiting solution for A in the form 

U A = Ei/i (r) + E, A (r) + E3/3 (r) + E,/, (r) + E j; {r,p) + EJ, {r, p)* . (76). 

^- I have not found it possible to investigate the form of these two last solutions in the general case, but 
it is easy to examine the nature of the solutions at infinity, when the nebula extends to infinity, and this 
enables us to form some idea as to the general nature of the solutions. Suppose that at infinity we have 

Ap dr dr 






in which a is real, then it can be shown that A - c^ (r, p) A^ &c., in which A', B', C, are functions of 
r only, and 

<i> (r, p) -= E5 cos (2 slan(n -^^T) r~'i^/isp) + Ee sin (2 ^/«m7^^T^) r"^/^/^?^) 

when the negative sign is taken in the above ambiguity, the circular functions being replaced by hyperbolic 
functions when the positive sign is taken. The value of ip is wholly real when squares of e> may be 
neglected (cf, § 13). 
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If we deduce the values of B and C from the solution (76), and substitute in the six 
boundary equations the values so obtained, we shall be left with six linear and 
homogeneous equations between the six E^s. Eliminating the six E^s, we have a 
single relation between ?^, the constants of the nebula and p. Now it will be seen 
that it will always be possible to pass to the limit ^ = in this equation, since this 
amounts only to finding the ratio of the values of /^ ^^ f% ^"t the two boundaries. 
The equation obtained in this manner will give us a knowledge of the configurations 
at which a change from stability to instability can take place. 

§ 19. It therefore appears that it is not sufficient to consider vibrations of frequency 
^ = as represented by positions of "limiting equilibrium." The method of 
PoiNCAR]^^ for determining points of transition from stability to instability is not 
sufficiently powerful for the present problem ; indeed it appears that it is liable to 
break down whenever there are boundary-equations to be satisfied, t 

It is of interest to notice that this complication is not (as might at first sight be 
suspected) a consequence of our having taken thermal conductivity into account. 
For we can put C = and remove the equation of conduction of heat without causing 
any change in our argument, except that the right-hand member in equation (74) must 
be replaced by a determinant consisting only of the minor of the bottom right-hand 
member in the present determinant. The value of A is now 

A = |)^ D^ + Dg, 

and the number of boundary-equations is of course »educed from six to four. . Thus 
an exactly similar situation presents itself, although we are now dealing with 
a strictly conservative system. 

The consequences of this result are more wide-reaching than would appear from 
the present problem, inasmuch as all problems of finding adjacent configurations of 
equilibrium are affected. For instance, it appears that an equilibrium theory of tides 
is meaningless except in very special cases (e.^., when the elements of the fluid in 
which the tide is raised are physically indistinguishable). 

If we attempt to calculate by the ordinary methods the tide raised in a mass of compressible fluid by 
a small tide-generating potential, we reach a number of equations which are (except in special cases) 
contradictory. To take a simple case, suppose we have a planet of radius Eq covered by an ocean of 
radius Ei, the whole being surrounded by an atmosphere which maintains a constant pressure w at the 
surface of the ocean. Let the law of compressibility be 'zr = c/o, where c varies from layer to layer of the 
ocean. Let the tide generating potential be %r^^Stt. Then the equations of this paper will hold if we 
write J? = 0, C = 0, ignore the equation of conduction of heat, replace XT everywhere by c, and include in 

^ " Sur TEquilibre d'une Masse fluide . . . ."— ' Acta. Math.,' 7, p. 259. 

t There is not, of course, a flaw in PoiNCARi's analysis, but he works on the supposition that the 
potential-function is a holomorphic function of the principal co-ordinates, and this supposition excludes a 
case like the present one. 

D 2 
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V a term aor*\ Equation (39) gives (except in the special case oi c = constant), B = 0. Equations (50) 
and (51) remain unaltered, and give a solution of the form 

A = Ei/i (r) + E2/2 (r). 

Now we must have A = when r = Ro, and this determines the ratio E1/E2. Also equation (49) must be 

satisfied, and this leads to a second and different value for E1/E2. 

A second example, of less interest but greater simplicity, will perhaps help to elucidate the matter. 

Imagine a non-gravitating medium in equilibrium under no forces inside a rigid boundary. Let the law 

connecting pressure and density for any particle be 'zr = Kp, where k varies from particle to particle. In 

equilibrium ^xr has a constant value wtq. Suppose now that we attempt to find an adjacent configuration 

which is one of equilibrium under a small disturbing potential V. The general equations of equilibrium 

are three of the form 

dV ^ 1(M 

dx p dx ' 

If the position of equilibrium only varies slightly from the initial position, dw/dx will be a small quantity 
of the* first order, so that (to the first order of small quantities) p may be replaced by its equilibrium value 



tro/K. We now have 



dw 'zro dV 



iX/X K Q/X 

and therefore, since ir is a single -valued function of position, 

1 ^V 



K ds 



ds = (i.), 



the integral being taken along any closed path. Since V and k are absolutely at our disposal, this 
equation is, in general, self contradictory. What we have proved is that there will only be an " adjacent " 
configuration of equilibrium under a potential V if V is a single valued function of k, a condition which 
will not in general be satisfied by arbitrary values of V and k. 

It is not difficult to see the physical interpretation of this last result. There were initially an infinite 
number of equilibrium positions, and therefore an infinite number of vibrations of frequency j9 == 0. To 
arrive at the configuration of equilibrium under the disturbing force we must imagine vibrations of 
frequency ^ = to take place until equation (i.) is satisfied ; the disturbed configuration will then differ 
only slightly from the configuration of equilibrium. For instance, if the disturbing field of force consists 
of a small vertical force g, the fluid must be supposed to arrange itself in horizontal layers of equal density, 
before we attempt to find the disturbed configuration. 

The interpretation of the result obtained in the first instance is similar, but more difficult. Consider 
a linear series of equilibrium configurations, obtained by the variation of some parameter a, such that the 
spherical configuration of our example is given by a = 0. The other configurations are not symmetrical, 
the asymmetry being maintained, if necessary, by an external field of force. Every degree of freedom in 
the configuration a = must have its counterpart in the configurations in which a is difterent from zero. 
In particular, the principal vibrations of § 12, in which (for the configuration a = 0) the dilatation, normal 
displacement, and temperature-increase all vanish, must have counterparts for all values of a. But when 
a is different from zero, the above three quantities cannot be supposed to all vanish. In general, therefore, 
these degrees of freedom provide solutions of the volume-equations, and these solutions contribute to the 
boundary-equations. In the special case of a = 0, these solutions do not affect the boandary-equations at 
all, so that to rectify the boandary-equations we must, so to speak, take an infinite amount of these 
solutions. In other words, the complete vibration of frequency p = becomes identical with one of the 
vibrations of § 12, in which u, A, and Ti all vanish. 
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An Isothermal Nebula. 

§ 20. Let us now examine the form assumed by our equations in the simple case in 
which X and T are the same at all points of the nebula. We find that, considering 
only the equations for the case ofp:=:: 0, equation (39) reduces to 

P — (7^\ 

and, in virtue of this simplification, the equation of conduction of heat (36), and the 
two thermal boundary conditions (33 and 34) are satisfied identically. We are left 
with equation (55) to be satisfied throughout the gas, and equations (32), (35), (52), 
and (53) to be satisfied at the boundaries. 

The solution of equation (55) is given in equation (56). Now we must satisfy 
equation (32) by taking B = at r = R^, and this, by equation (50), gives the value 
of A at r = ~Ri in terms of Ej and E^. Hence equation (52) reduces to a homo- 
geneous linear equation between E^ and E^. 

When n is different from unity, we satisfy equation (35) by taking A = at 
r = Rq, and this reduces equation (53) to a homogeneous linear equation between 

joj-i and JlI^o. 

When n = 1, equation (35) reduces to a linear equation between (A)^^^^, E^ and E^. 
Equation (53) is a second equation of the same form, and the elimination of (A),,^r^ 
from these two equations leads to a homogeneous linear equation between E^ and Eg. 

Thus, in either case, we see that the whole system of equations reduces to a pair 
of homogeneous linear equations between Ej^ and E^. The elimination of these 
quantities leaves us with a single equation between n and the constants of the 
nebula. 

We can, therefore, satisfy all the equations for a vibration of frequency p -= by 
imposing a single relation upon the constants of the nebula. The unknown solutions 
which are multipHed by Eg and Eg have not been taken into account at all, but since 
the condition that there shall be a vibration of frequency p = must of necessity 
reduce to a single equation, it will be clear that if these solutions had been taken 
into account, we should have found it necessary to take E^ = Eg = 0. 

Thus, in the case which we are now considering, a vibration of frequency p = 
is equivalent to a configuration of limiting equilibrium. It is not hard to see that 
this results from the fact that the particles of which the nebula is composed are 
physically indistinguishable. This very fact, however, introduces a further complica- 
tion into the question. It will be noticed that, although the value of f has been 
found at every point of the nebula, it is impossible to determine the separate values 
of A and B. On the other hand, the physical vibration must have a definite limiting 
form when p = 0. Now it is easy to see that a motion of the gas in which 
f vanishes at every point of the gas, and in which A and B vanish separately at the 
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boundary, will, in every configuration of the gas, satisfy our equations with p = 0, 
Such a motion, in fact, simply leads to a configuration which is physically 
indistinguishable from the initial configuration, and in which the potential energy 
remains unaltered. The motion which we have found from our equations is the sum 
of a motion of this kind, and a true limiting vibration. It is impossible to separate 
the two motions, except by considering vibrations of frequency different from zero, 
but fortunately the question is not one of any importance. 

§ 21, Let us now attempt to form the final equation in some cases of interest. The 
equations of an isothermal nebula at rest under its own gravitation have been 
discussed by Professor Darwin.^ Our function te {equation 54) is given, in the case 
in which the nebula is isothermal, by the equation 

and it will be seen that this is the same as the ii of Professor Darwin's paper. It 
appears that in general u cannot be expressed as a function of r in finite terms, but a 
table of numerical values of u is given, f The value of ii approximates asymptotically 
to unity at infinity, so that at infinity p varies as ? ~^ Darwin^s nebula extends 
from r = to r =: co ^ bu.t it is obvious that we may, v/ithout disturbing the 
equilibrium, replace that part of the nebula which extends from r = to r ™ Rq by 
a solid core of mass equal to that of the gas v/hich it replaces. We may also remove 
that part of die nebu.la which extends from t' :rz R^ to r === co , if we suppose a pressure 
to act upon the surface r = R^ of amount equal to the pressure of the gas at this 
surface. We may suppose the medium outside this surface to be of any kind we 
please, but as it has already been pointed out that the pressure can, in nature, only be 
maintained by the impact of matter, we shall suppose that this matter is of a density a 
which is continuous with the density p of the nebula at the sin^face of separation. 
We may now write equation (52) in the simple form 



i i^'''-) 



d) 



f=:R, 



(7^) 

\J O h 9 9 9 A 9 9 9 \S't//t 



We have, up to the present, supposed the nebula to be acted upon by a spherically 
symmetrical system of forces in addition to its own gravitation. Now it is essential 
to the plan of our investigation that we shall be able to make the configuration of the 
nebula vary in some continuous manner, and this compels us to retain this generali- 
sation. We shall, however, suppose that when the nebula extends to infinity, 
It retains some definite limiting value t^^, thus including the free nebula as a 
special case. 

* G..H. Darwin, 'Pliil. Trans.,' A, vol. 180, p. L 
t Loc, cif,, p. 15. 
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§ 22. Let us, in the first place, consider the "series" of nebuliB such that u has a 
different constant value for each. This series includes a single free nebula, for it 
appears from Darwin's paper that there is a nebula such that '?^ = 1 at every point. 
This nebula, it is true, has infinite density at the centre, but this objection disappears 
when the innermost shells of gas are replaced by a solid core, the mean density of the 
core being equal to three times the density of the gas at its surface, and therefore 
finite. Let us, in the first instance, simplify the problem by supposing that the core 
is held at rest in space. The boundary equations (35 and 53) which have to be 
satisfied at r = Rq now take the forms 

\A./r=Ro ^^^ ^ ......... (SOj, 

-~y^^{gr "" ) = (81), 

independently of the value of n. The value of u in equation (55) being now 
independent of r, we may write the solution (56) in the form 

in which /^, fi' are the roots of the quadratic, 

^ (^ — 1) — f^ (ti + 1) •>- 2ti . . . . . . . . (83). 

We accordingly have 

IJL -{fju' :=z I ; ^ — jLt' =: 2^/(n + 2")^ "^^ ^'^ 5 i^/^' = "" ^^ ('^^ + 1) + 2u ., (84). 

Equation (79) now takes the form 



El (/x + n) -Rf^--^ + Ec, {yJ + n) Rf '-^^™^ =:; .... (85), 
while equation (81) becomes 

E^ (^ -^ n ^ 1) Ro^-"'^"^ + E^ (ft' -n-l) Ro^'-^^"^ = . . . (8G). 
The elimination of E| and E^ from these equations gives 

Rq/ {fjb + n) {^jb — n — 1) ' * * * ' ' * ^ ^* 

The fraction on the right hand can be simplified by the help of equations (84); it 
is equal to 

^ (' ^ — 0^ + if) + (/^ — ^') (^^ + i ) 
2'(u ~(n+ if) ^ (/. - ^) (^ + i) * 

Now the left-hand member of (87) may be replaced by 

cosh {| {fjb — fjb^) log (E^/R())} + siiih {|- (/i ~ ^0 log (Ki/^o)} 
cosh { H/^ - /^O log (Ki/Ko)} - sinh {|(/i - //) log (Ri/Eo)} 
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so that the equation itself reduces to 

tanh {"I (fi — /x^) log (R|/Ro)i 



(f6 — /^') (11 f D 



w — (ri 4- |-) 



» e « 



(88). 



This equation expresses the relation which must exist between II|/IIq and fi — fi^ 
(or, what is the same thing, between Ri/Rq ^^^^^ '^^)j i^^ order that p === may be a 
solution of the frequency equation. 

§ 23. We shall be able to interpret this equation most easily by adopting a 
graphical treatment. If we write 



X 



4 l/^ /x j , y^-^ 



X -{• in ~\- ^y 






Ih 



jj 



tanh { ^x log (Rj/Rq)}, 



then the equation can be written in the form 



Vi = y%^ 



It will be noticed that y^ remains real when x is negative, an equivalent expression 

for 1/2 being 

1 „__„..^ 

y^ = -^^y=r^ tan { y/ - a; log (Ri/Rq)}. 



(A/ 



The roots of equation (87) are now represented by the intersections of the graphs 
which are obtained by plotting out yi and ?/g as functions of x. These two graphs 




y 



o 



Fig. 1. 



i-C' 




Fig. 2. 



are given in figs. I and 2 respectively, the graphs being drawn separately for the 
gake of clearness. The graph for y^ is, of course, the same for all values of R^/R^ ; 
that for y^ can be varied so as to suit any value of Ri/Rq by supposing it subjected 
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to ail appropriate uniform extension parallel to the axis of y, and contraction parallel 
to the axis of x^ or vice versa. Similarly, different values of (^^ + g) can be 
represented by contraction and extension of the first graph. 

If we imagine these two graphs superposed, we see that there cannot, under any 
circumstances, be an intersection in the region in which x is positive, ^.e. (equation 
(84)), for a value of u less than -| (n + ^y. The lowest value of u for which an inter- 
section can possibly occur is u = 1, and this occurs only when Ei/Rq = oo . As Rj/Ro 
decreases from infinity downwards, the lowest value of tt for which an intersection 
occurs will continually increase. Whatever the value of R^/Rq may be, there are 
always an infinite number of intersections in the region in which t^ > ^ (n + ^y. 

The values of u found in this way determine the ^^ points of bifurcation" on the 
linear series obtained by causing ti to vary continuously. Thus we have seen that as 
It continually increases the first point of bifurcation of order n is reached when it has 
a value which is always greater than ^ (n + -kY- When Rj/Rq is very large, the first 
point of bifurcation is of order ^ = 1, and its position is given by 

7/ -— 1 1 (Qq) 

W " '■ J. Q» « s • t » . « • a « \ \J %J J i 

§ 24. Let us, in future, confine our attention to the case in which Ri/Rq is very 
large. If we gradually remove the restriction that u is to be independent of r, the 
various vibrations of frequency p = will vary in a continuous manner. Equation 
(55) remains unaltered in form, and, at infinity, it assumes the definite limiting form 

r"^ A = {n(7i + 1) — 2i^^J'^^ , (90), 

where u^ is the limit (supposed definite) of it at infinity. It therefore appears that 
at infinity the solution for ^ approximates asymptotically to that given by equation 
(82), if /x, fji^ are now taken to be the roots of 

^^ — . 1) = n{7i + 1) — 2u^ ...*... (91), 

Equation (85) accordingly remains unaltered. Equation (81) takes a form which 
is no longer represented by equation (86), but which will impose some definite ratio 
upon Ej and Eo. It is therefore clear that when Rj^ is very great, equation (85) can 
only be satisfied, at any rate so long as /x and /x' are real, by taking /x — ja^ very 
small. Thus a point of bifurcation will again be given by /x — /x^ = 0, our previous 
investigation sufiicing to show that this gives a genuine solution to our problem, and 
does not correspond to an irrelevant factor introduced in the transformation of our 
equations. This point of bifurcation is moreover the first one reached as u increases, 
since it is at the point at which /x, /x^ change from being real to being complex. 
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We conclude that, independently of tlie values of lo at points inside the nebula, 
the smallest value of u^ for which a vibration of zero frequency and of order n is 
possible is given by 



'«■- ^ i {n + ^y 



(92). 



or, for all orders, is given by 



II 



OC' 



J c 



8 



(93). 



tlie limitiiig vibration being of order n ^ J. 

It ought to be noticed that for this limiting Adbration equation (82) tails to 
represent the solution owing to jjl and jjl^ becoming identical. The true solutions for 
real, zero, and imaginary values of /x — /.t' may be put respectively in the forms 

f = Cj \/r sinh {^{(jl — /x^) log eR}, 



^=1= CjVrhgeR, 



f = Civ/"/' sin < ^ (/x — /x^) log ell 



in which Cj^ and e are constants of integration. 

At infinity p vanishes to the order of 1/r^, so that dp/dr 
^ for very great values of r is therefore (equation (50)) 

^=:=XT(-»' 2A+ B4 



2p/'r\ The value of 



At the outer boinidary a surface ~ equation (32) directs us to take B -^ 0. 
Following this out, we fi.nd that at infinity A is of the same order as ^, and therefore 
becomes infinite to the order of ^/r. Suppose^ on the other hand, that we start by 
taking A = 0, so that B = f/XTr. The value of B now vanishes at infinity to the 
order of l/\/r, and the surface-equation (32) is satisfied by a motion which vanishes 
at infinity. It would therefore appear to be easier to satisfy the boundary conditions 
when r is actually infinite than when r is merely very great. This result opens up a 
somewhat difficult question, which will be considered in the next section. 

Before passing on, we may consider in what way the results which liave already 
been obtained will be modified, if we suppose the core of the nebula to l^e free to 
move in space, instead of being held fast. For the free nebula ii^ = 1, so that our 
results show that a free nebula will be stable if the core is supposed fixed. The 
same must therefore obviously be true when the core is free to move, since a motion 
in which nebula and core move as a single rigid body will not influence the potential 
energy. When the nebula is not free, fixing the core may be regarded as imposing 
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a constraint which does no work ; freedom of the core therefore tends towards 
instabihty. It will be proved in § 28, that a nebula is stable for values of ii^ which 
are less than the critical value, and unstable for values greater than this value. 
Assuming this for the moment, we see that a nebula in which the core is free to move 
must necessarily be unstable if n^ has a value greater than 1 1. 

If then, we start with a free nebula and imagine n^ to gradually increase from 
11^ ^ 1 upwards, the core being free^ it follows that the nebula will first become 
unstable when n^ reaches some value sucli that 



Ig- j> u^ ^ I 



. . , (94). 



§ 25, Tlie nebula extending to infinity, let us attempt to find the displacement 
which will be caused by a small disturbing potential v,^ given by 



Vn 



inr 



CL 







2n + 1 Ir^^ + i ^ 






(95). 



It is clear that the displacement required will be given by our equations if we 
include in V (equation (27)) the terms 



4 



TT 



2??/ + 1 






The equation replacing (42) may be transformed in the manner of § 15, and the 
resulting equations will be those of § 16, except that we must replace (52) by 



1 d 






Ch 



a, 



"A (p - a J 



i"=R 



v, 



r 



n — l 



:R, 



(96), 



and (53) by a similar equation. 

If a displacement can be found to satisfy these modified equations, the external 
disturbing potential which will be required to hold the system in this dis23laced 
position will be given by equation (94). Now the condition that this displaced 
position shall be one of limiting equilibrium is that this disturbing potential must 
vanish. To be more precise, v,i must be such that the force derived from it vanishes 
at every point of the nebula. We must therefore have 



ai7 



,n 







at all points of the nebula, including r = R^ Now (95) may be regarded as an 
equation giving a^ in terms of R^. Taking p =orj, as before, we find from (95) 



^1''*' 



r 



H 



' 1 d , ^ J 
R « dr ^^ ^ 



r=K. 



'u 9 
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and this vanishes at all j)oints, including r = Rj, in the case in which Rj is put equal 
to infinity^ if 

■X~*U •— «» \j ^ ^ ^ ^ ^ ^ ^ ^ J J \ tJ f J t 

I 



r =<xi 



The condition that ciqIv"''^^ shall vanish at every point would lead to a similar 
equation to be satisfied at the origin, if there were no core. If, however, we retain 
the core, it leads to the same equation as was found in § 22 (equation (81), when the 
core is held at rest). Thus, our present method of finding a position of limiting 
equilibrium has led to a result different from that ol)tained by the search for a 
vibration of zero frequency^ in that equation (97) replaces equation (79). 

The value of f at infinity is given by equation (82) ; hence w^e have 

U ^- = rE.r'^-"^ + 'E,r^''"'\^^ . \ . . / \ . . (98). 

;' :^ 00 ' 

As before, the equation to be satistied at r r-z E^^ determines tlie rtitio of E| to E,> : 
equation (97) is therefore satisfied, if the i*eal parts of /x and. /x^ ai'e each less than 
unity. Now jul, fju^ are the roots of equation (91), hence this condition is satisfied 
provided 

n(n + 1) < 2t^^ . . . . . . . . . (99). 

§ 26. Let the kinetic and potential energies of a small displacement be given^ in 
terms of the principal co-ordinates, by 

2 V == ^x^i + &2'%^ + . . » 



so that the equations of motion are 



a|X| — &;}_x = 



&c., and p^ is given by 

The method of §§ 20-24 amounted to finding vibrations such, that j^'^ = 0, and 
therefore^ by equation (100)^ solutions of 

In § 25, on the other hand^ we started with the supposition that the nebula extended 
to infinity, so that all the quantities a and 6 are liable to become infinite. The 
equation giving vibrations of frequency "p. ~- is no longer equation (101), but is 

T « _1. — A /I AOX 

and thjs is obviorisly .more general than equation (101), 
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It will be noticed that the method of §§ 20-24 is the method which is mathematically 
appropriate to the .case of a nebula enclosed in a surface maintained at constant 
pressure, while the method of § 25 is that appropriate to an infinite nebula. In the 
former case, a vibration of frequency p =: may represent a real change from stability 
to instability ; in the latter case such a vibration leads to an adjacent configuration of 
equilibrium, and is, in this sense, a point of bifurcation, but does not denote a change 
in the sign of j/. 



The General Case of a Nebula extending to Infinity, 

§ 27. The method to be followed has been explained in § 18. The general 
diffei'ential equation is of the sixth order. Four solutions have definite limiting forms 
when j9 = ; the remainbig two take singular forms. The former have been examined 
in I 16 ; tlie latter are represented mathematically (p. 18) by functions which do not 
approach a definite limit as p approaches a zero value, and physically (p. 11) by 
systems of steady currents. 

There are six constants of integration, E^, E^, E3, E4, E^, Eg, of which the two last 
belong to the singular solutions. Let us suppose (as is always possible (p. 19)) that 
the ratios of these six constants are determined from five of the boundary-equations, 
that which is not used being the equation satisfied by ^ at the outer boundary. This 
remaining boundary-equation now takes the form {cj\ equation (56)) 

Eit/.i(R,) + Eot/.,(EO + E,t/;,(EO + Eet/;o(R^^^ . . . (103), 

in which the four E's are definite quantities. The four ^\fB must have definite 
limiting values (zero and infinity being included as possible values) when R^ = 00 . 
Thus in equation (103) some terms must preponderate over the others. When the 
nebula is isothermal, these terms are the first two. Hence, when the nebula is not 
isothermal, it follows from the principle of continuity, that the same two terms nuist 
still preponderate, at any rate for some finite domain including the isothermal 
nebula. Otherwise it would be possible to change the stability or instability of 
a nebula by an infinitesimal change in the physical constitution of the nebula. 
Hence throughout this domain, equation (103) must reduce to its first two terms, 
i.e,, must become formally the same as in the case of the isothermal nebula. But 
the solution for f (and therefore the functions t//^, i//^), remain formally the same in the 
general case as in this particular case, and therefore the stability-criterion derived 
from equation (103) remains formally the same. 

It follows that whether the nebula is isothermal or not (provided always that the 
configuration lies within a certain domain of equilibrium configurations) the critical 
configurations are given by the two equations (92) and (99). 



1 



aO M\L J, IL JEANS ON Till: STABILITY OF A HPHElilCAL NEBULA. 



Exchange of Stahilities, 

§ 28. We have now completed an investigation of the conligurations at which 
a transition from stabihty to instability can occur, as regards the spherical form, for 
vibrations of orders clifierent from zero. It is unnecessary to discuss vibrations of 
order n ^^ 0^ for the following reason. 

Our problem is to determine the changes in the configuration of a nebula which 
will take place as the nebula cools^ startnig from a spherical configuration, supposed 
stable. We are not concerned with the succession of spherical configurations, but 
only with an investigation of the conditions under Avhich a spherical configuration 
becomes a physical impossibility. Now a point of bifm^cation of order n ^ does 
not indicate a departure from the spherical configuration. It indicates a choice of 
two paths, one stable and the other unstable, and the configurations on both paths 
will remain spherically symmetrical. 

We have therefore determined already the circumstances under which a transition 
from a symmetrical to an unsymmetrical configuration can occur. It remains to show 
that there is, in effect, an exchange of stabilities at a point of bifurcation, and to 
examine on which side of the point of bifurcation the spherical configuration is stable. 

We are going to prove that the spherical configuration is stable for all values of u 
less than Uq, the lowest value of tt at which a point of bifurcation of order different 
from zero can occur. Our method will be as follows : Any two equilibrium configura- 
tions can be connected by a continuous linear series of equilibrium, configurations, and 
ti will vary continuously as we move along this series. If one of the two terminal 
configurations is stable, and if the linear series can be chosen so that ?i does not at 
any point of it pass thiTjugh a value for which a vibration, of frequency p ™ is 
possible, then we know that the other terminal configuration is also stable. 

The value of y, the gravitation constant, has been taken equal to unity. If this 
constant is restored, the value of ti becomes (equation (54)) 



Vi 



o dp I dm 

' ' dr / dr 



ince 'd '==^ XTp, we have 



dm d /s^rr\\ I \ m dp 

- rr:: n -- I Ai ) + Al " 

dr ' dr ^ ^ dr 



For an infinite nebula, the first term on the rightdiand of this equation will 
vanisli at infinity in comparison with the seconcl Hence we have as the 
value of tt^ 
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If we write y = we pass to the case of a non-gravitating nebula, and we see 
that ?/ ^ = provided the ratio of pr^ to XT remains finite at infinity. Now we can 
keep the vahie of p and XT the same at every point by subjecting the nebula to 
an appropriate external field of force, and this field of force will be exactly the same 
as the gravitational field which was annihilated upon putting y =: 0. It is spheri- 
cally symmetrical, and its potential vanishes at infinity to the order of 1/r, so that 
it comes within the scope of our previous analysis. For values of y intermediate 
between the natural value (y =: l) and the value y =:^ we can obtain the same 
result by taking a field of force equal to 1 — y times the foregoing. As we increase 
y from to 1 we obtain a linear series, in which the configuration of the nebula 
is unaltered, the nebula being gradually endowed Avith the power of gravitation. 

For the general configuration of this series, consider the work done in a specified 
displacement, which is proportional to S^^, at every point. The potential (gravitational 
-f that of external field) after displacement Avill be of the form 

Avhere a and h are functions of r and independent of y. The total work done against 
this field during the displacement is therefore of the form 

By, 

where B is independent of y and depends solely upon the particular displacement 
selected. The work done against the elastic forces is of course independent of y, 
and depends solely upon the displacement selected. This work is essentially positive. 
The total work is therefore of the form 

A + By, 

where A is positive and B may (§ 2) be negative. Since y is proportional to u^ this 

may be written 

A + B'?i^^ (105). 

Suppose this function calculated for all possible displacements. Then we shall find 
that for values greater than some definite value of u^ it is possible for the work 
done to become negative. For values of u^ less than, this critical value, the work 
will be positive for all displacements. Hence from the form of expression (105) 
it follows that the passage of ii^ through a critical value denotes a real change from 
stability to instability, and that the stable configurations are given by the smaller 
values of u^^. 

Meccqntulation and Discussion of Results. 

§ 29. We have seen that the vibrations of any spherical nebula may be classified 
into vibrations of orders n z= 0, 1, 2, &c., a vibration of any order n being such that 
the displacement and change in temperature at any point are each proportional to 
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some spherical surface harmonic S,, of order n. The frequency of vibration is 
independent of the particular spherical harmonic chosen, depending only upon the 
order n. 

The vibrations of order n = have been seen to be of no importance ; the stability 
of the vibrations of orders different from zero has been discussed, in the limiting case 
in which the nebula extends to infinity, with the following results :— 

Starting from any stable configuration of spherical symmetry, the vibrations of any 
order ri, different from zero, all remain stable until the function u^^ defined by 
equation (104), passes through a certain critical value. In any case this critical 
value is first attained for a vibration of order n = 1. 

For a nebula which actually extends to infinity, the critical value is ii^^ ^^ L 
When this value is reached we come to a second series of equilibrium configurations, 
tlie form of which will be investigated later. If this value is passed, the configura- 
tion remaining spherical, there will not be vibrations in which the time enters 
through a real exponential factor, because the critical vibrations remain of frequency 
p ::=:= 0, the inertia of the nebula being infinite. 

If the radius 11^ of the nebula is regarded as very great but not infinite^ this 
statement is not true, since the inertia cannot now become infinite. In this case the 
first new series of equilibrium configurations is again reached when {ii)^^^ji>^ attains a 
certain critical value, and the ciitical vibration is again of order n =^ 1. The critical 
value of {u)y^-^^ has not been calculated, but when 11^ becomes infinite, it has a 
limiting value which has been shown to lie between ] and 1|-. 

Taking y =: 1, we have as the value of tt,,, 






■-/'::= 00 



The question of stability turns entirely upon the value of this function, which may 
appropriately be termed the ^' stability-function^" 

We now see that the whole question of stability depends upon the ratio of the 
density to the elasticity at infinity. This result is not hard to understand. In the 
first place, since the nebula extends to infinity, we may^ so to speak, measure it upon 
any linear scale we like. If we measure it on a sufficiently great scale, tlie nebula 
still remains of infinite extent, but tlie variations in temperature o:i' structure which 
occur near the centre can be made to appear as smali as we wish, and the solid core 
can be made to appear as insignificant as we wish. Thus by measuring any nebula 
upon a sufficiently great scale we can make it appear indistinguishable from, an 
isothermal nebula, and the critical vibration for which p == does not disappear from, 
sight, since in the limit this vibration (measured by f/r) remains finite at infinity. 
Further, as Professor Daewin points out, we can make it appear like a nebula in 
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which u maintains a constant value throughout.^ Passing on, we notice that the 
stabihty function now depends solely upon the ratio of the density to the elasticity. 
The different elements of the nebula are attracted towards one another by their 
mutual gravitation, and are kept apart by the elasticity of the gas. For certain 
values of the ratio of these two systems of forces, it will be possible to find displace- 
ments in which the work done by one system exactly balances that done against the 
other, and these are the critical vibrations. 

The stability function it^ is a function only of the quantities determining the 
equilibrium configuration of the nebula, and its value may therefore be found from 
the equations of equilibrium. We proceed to examine the value. 

Evaluation of the Stability Function. 

General Case of a Nebula at Rest 

§ 30. We have already quoted Professor Darwin's result that 1/^ = 1 for an 
isothermal nebula at rest, and the considerations put forward in the last section will 
probably suggest that the result in the more ge^neral case will be found to be 
independent of variations in temperature at finite distances, provided only that 
the temperature has a definite limit at infinity. We shall, however, examine the 
question ah initio^ using a slight modification of Darwin's method, and making the 
problem more general by retaining a spherically symmetrical system of external 
forces. 

We shall denote the potential of this system of forces by V, and use V to denote 
the gravitational potential of the nebula itself The total potential is now V + ^^ 
so that the equation of equilibrium, equation (11), takes the form 

' '^ (xTp) - I: (V + V) = 0, 



p dr ^ ^^ dr 
and if M is the mass of the solid core, this can be written 



^"|(XTp) + 47r[^/.r^^r + M-r^^' = .... (107). 



Difierentiating with respect to r, 



" ' " " (XIri) + 4.p,- - i: (r' Z) = 0. 



dr \ p dr ^ ^' I ^ dr \ dr j 

Write 

XTp = e^, 

^ Gr. H. Darwin (^.c ante, p. 16), " If we view the nebula from a very great distance, . . . the solution 
of the problem becomes y = log 2a:2." Now u = ~ ^x^ d^y/dx'^, so that this solution is equivalent 
to u = 1. This justifies our statement, and shows at the same time that for any nebula at rest and 
in equilibrium u^ has the critical value u^ = 1, provided it is acted upon by no forces except its own 
gravitation. 
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and 

Too ^^. 

np ~ — I s 

SO that 

d^ _ 1 d 

dr ~~ 'xT>3 f^ig^^ ' 

then the above equation may be written 

dhf d I dY^\ 

y{ + 47reV - T^ ,- r^^V =0 (I08l 

rfa-" dr \ dr / ^ ' 

At infinity we are supposing XT to have a definite and finite limit, so that the 

limiting value of x is l/XTr. Let us further suppose that 4-~-( r^-^ ) has a definite 

-'■ ^ dj' \ dr I 

limit given by 

dr\ dr ^ ^■^' 

and that squares of V^^ may be neglected. Then tlie limiting form of (108) at 

infinity is 

dh/ ^irey Y'' 

Write 
then 

dhi d?-r) 2 d? , , ,^, 

and 



Equation (110) is now transformed into 

S + i(^^~l) + 4|,log(XT)~-^^^^, = . , . . (111). 

d^ 
In the special case in which , ^^ log (XT) vanishes, this may be written 



d',. 



x^ 



and at infinity (^.e., for very small values of ^), the solution is 



77 



+ /\/4--cos(|-y7 log^ j , . . . . . (112) 



where A, B are the two constants of integration, 
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In the more general case in which — log (XT) cannot be supposed to vanish, it is 

LttAJ 

clear that this term will vanish at infinity in comparison with the other terms 
in (HI), if 7^ has the limiting value given by (112), and therefore that (112) is the 
limit, at infinity, of the solution of (108). 

Of the two arbitrary constants, A and B, the former corresponds to the indeter- 
minateness of the linear scale upon which the nebula is measured, the second to the 
indeterminateness of the conditions at the inner surface of the nebula. If there is no 
core, there is only one value of A/B which will give a finite density of matter at the 
centre of the nebula. Further information as to equilibrium configurations can be 
found in Professor Darwin's paper, or in a paper by A. Bitter. "^^^ 

For our purpose it is sufiicient to know that the second term in t) vanishes with x 
for all values of A and B. Hence at infinity 

t/ = -4- 00" 



■\3T3..<3 X'T 



XT 27!- 277?'^ 



5 



and hence (equation (107)) 



Putting V'^ =: 0, we arrive at the anticipated result that the stability function has 
a unit value, for every nebula which extends to infinity in such a way that XT 
has a finite limit at infinity. 



A Sloivly Rotating Nebula, 

§ 31. The case which is of the greatest physical interest, is that in which the 
nebula is not at rest but is rotating in a position of relative equilibrium 

Here the arrangement is no longer in spherical shells, so that the foregoing analysis 
breaks down. If, however, we suppose the rotation o) to be so small that o)^ may be 

J. a- J ■ 

neglected, it will be easy to modify the foregoing analysis, so as to take account of 
rotation. 

We shall still suppose the nebula to extend to infinity, so that v^^e must not suppose 
the rotation to be the same at all distances, for in this case a finite value of ca would 
imply an infinite velocity of those parts of the nebula which are at infinity. Let us 

* ' Wied. Ann./ vol. 16, p. 166, 

F 2 
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suppose that at infinity the Hnear velocity approximates to a finite limit, so that we 
may write 

for all values of r greater than a certain amount.* 

So long as we are only concerned with configurations of equilibrium and vibrations 
of frequency p == 0, the rotation may be allowed for by the introduction of a force of 
amount w V sin 6 per unit mass, acting perpendicular to the axis of rotation ; or, what 
comes to the same thing, by the introduction of a potential 



2 
'3 



. (1 — P^)^ oi%'c/r 



or 

where, for all values of r greater than a certain value. 

Let us examine separately the two effects arising from the two terms of this 
potential, beginning with the term — P^V, There will in this case be a correction 
to be applied to all equations, and this correction will consist of the addition of 
a small term containing co-^Pg. Let us suppose that all symbols which have so far 
denoted functions of r, denote in future the mean value of the corresponding 
quantities averaged over a sphere of radius n For instance, p is no longer the 
density at distance r from the centre, but is the mean density over the sphere of 
radius r. The density at any point will be of the form p + ^^P^Pa? where p^ is 
a function of r. We may in every case equate the coefficients of different harmonics, 
and by equating the coefficients of terms which do not contain the terms w^P^, we 
shall obtain the same equations as were obtained in the case of w = 0, except that 
the meaning of every term is altered. 

The equations derived from the parts which do not contain co vnll suffice, as before, 
to determine p, so that the values of p are of the same form as before, except that 
the quantities involved have a slightly different meaning. Hence the stability 
criterion is still given by the value of the stability function u^ ; while equation (107) 

^ This particular law is chosen for examination because it leads most quickly to the required result. 
The case in which w vanishes at infinity more rapidly than l/r is covered by taking = 0. Heie, 
however, the angular momentum vanishes in comparison with the mass, and it is not surprising to find 
that a rotation of this kind does not aifect the question of stabitity. The case in which w vanishes less 
rapidly than l/r is physically impossible, since it gives an infinite linear velocity at infinity, but may be 
theoretically included in the case of = go . 

Any special assumption about the value of w at infinity must, however, disappear when we turn to the 
case of a finite nebula (§ 26), in which may be appropriately supposed to correspond to the surface 
velocity wKi. 
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remains true, if the new meaning is given to the symbols in each case. We conclude 
that the question of stability is not affected by the potential — Ps^'. 

The remaining potential term is the spherically symmetrical term Y\ The total 
potential may now be taken to be V+V', and this potential, besides being spherically 
symmetrical, satisfies the condition which was postulated in the determination of the 
criterion of stability ; namely, that its radial differential coefiicient shall vanish at 
infinity to the order of 1/r. The value of the derived function Y'^ (equation (109)) is 

Y" := U -~ ('1 2 '^) ^ |-^2^ by equation (] 1 4). 

Hence the stability function is given by (of. equation (113)) 

. ^^ 



3X3T3 • 



We have therefore found that when an infinite nebula is rotating, with such 
angular velocities that the linear velocities at infinity have the limiting value XI, the 
value of ic^ is greater than unity no matter how small H may be. This result has 
only been obtained on the supposition that co^ may be neglected. We have obtained 
no information as to what happens when co^ is taken into account, ^.e., when the 
square of the '^ ellipticity " of the nebula is taken into account. 

Influence of Viscosity, 

§ 32. No account has so far been taken of the viscosity of the gas. The terms 
arising from viscosity which may be supposed to occur in the true equations of 
motion, will contain the coejfficient of viscosity (fi), and will in each case- depend on 
velocities and not on displacements. Hence viscosity enters the equations of motion 
through the factor [jbip. The vibrations for which ^ = are accordingly unaffected 
by viscosity, and since it is upon the existence of such vibrations that the whole 
question of stability turns, it is clear that the results already obtained must remain 
true even in the presence of viscosity. 

It can be shown that equations (24) to (26) specify a principal vibration, whether 
the gas is viscous or not. The result is stated without proof, as the proof is rather 
lengthy, and has no bearing upon the main question under discussion. 

A Nebula in Process of Cooling, 

§ 33. In the mathematical investigation we have been concerned with vibrations 
about a position of absolute equilibrium. In nature, no such position of absolute 
equilibrium will occur ; the condition of the nebula will be incessantly changing. 
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Let us suppose the temperature of the nebula to be continually cooling, owing 
either to radiation of heat from its surface or to a process of quasi-evaporation such 
as is described in Professor Darwin's paper (§13 or p, 66), Since the gas (or quasi- 
gas) is not a perfect conductor, the nebula will not at any time be in perfect thermal 
equilibrium. The changes in density of all parts, and in the temperature of the 
inner parts of the nebula will, so to speak, lag behind their equilibrium values as 
determined by the changes in the temperature of the outer part of the nebula. It 
is, therefore, clear that so long as the nebula is cooling, the ratio of the density to 
elasticity in the outermost layers of gas will be greater than that calculated upon 
the assumption of perfect equilibrium. This ''lag" accordingly decreases the value 
of the stability- function, and so supplies a factor which tends to instability. 

Summary and Discussion of Results. 

§ 34. Let us now examine to what extent we have found solutions of the two 
problems propounded in § 4. 

Firstly, as regards the stability of a spherical nebula of very great size, of which 
the outer surface is maintained at constant pressure. We have found that the 
stability "function for such a nebula (in the limiting case in which the outer radius is 
infinite) has a unit value when the nebula is in equilibrium and at rest. This value 
is increased by allowing for the '' lag" in tempei'ature caused by the cooling of the 
nebula. It is also increased by a rotation of the nebula, at any rate so long as this 
rotation is small. The nebula will become unstable as soon as the stability-function 
becomes greater than a certain value, which has not been calculated, but is known to 
be between 1 and 1^. The investigation of § 23 leads u.s to expect that the critical 
value of the stability-function will increase as Rj decreases, although this has only 
been strictly proved for a single case. 

It is therefore possible that, even when the nebula is non-rotating, the temperature- 
lag may be sufficient to make the nebula unstable. If we disregard the temperature- 
lag, it seems probable that a small rotation will suffice to bring about instability. 
This latter question, however, deserves more detailed examination. 

§ 35. Let us suppose that the nebula starts from rest in a configuration of absolute 
equilibrium, and that the rotation is gradually increased. In this way we obtain a 
linear series of configurations of relative equilibrium. When the rotation is small, 
the configuration, instead of being strictly spherical is slightly spheroidal The 
series we are considering is therefore the analogue of the series of Maclaurik 
spheroids of an incompressible fluid. So long as the rotation remains small, we may 
separate the two terms of the rotation-potential in the manner explained in § 31. 
We may, in fact, suppose our analysis still to apply as if the configuration remained 
spherical, and the only efifect of the rotation is to increase the value of the stability- 
function. For larger values of co'^^ all our results are subject to a correction of the 
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order of a)'^\ For small values of co^, the value of co^ will be proportional as we have 
seen (§ 31) to u^-— 1, so that this correction may be supposed to be proportional to 
{u^—iy. The first points of bifurcation of orders 1, 2 occur (in the spherical 
configuration) at u^— 1 = 0, 2|- respectively, where 6 is known to be less than ^. 
Now it would seem to be fairly safe to neglect 6^, but even if we waive this point, it 
will be admitted that the correction of the order of {ii^ — 1 Y cannot be so great as to 
change the order in which these two points of bifurcation will occur. 

We therefore see that a rotating nebula will become unstable for a comparatively 
small value of oj^, the critical vibration being of order n ~ 1. The new linear series 
is one in which (except for the spheroidal deformation caused by the rotation) the 
surfaces of equal density remain spheres, which are no longer concentric. The linear 
series of order n =- 2 will accordingly be unstable : this is the analogue to the series 
of Jacobian ellipsoids in the incompressible fluid. 

§ 36. The case of a nebula which actually extends to infinity is much simpler. 
Here the value of u^ is again unity, and this value is increased, as before, either by 
fcemperature4ag or rotation. Every point at which u^ is greater than unity is in one 
sense a point of bifurcation, since starting from this point there is a series of 
unsymmetrical equilibrium configurations. Strictly speaking, these points do not 
indicate an exchange of stabilities, for the critical vibrations remain of frequency 
p = even after passing the point. They possess, however, the property that a 
critical vibration, if once started, will continue increasing, since the forces of 
restitution (of whichever sign) vanish in comparison with the momentum of the 
vibration. 

§ 37. Let us now try and examine which of these two hypotheses is best capable of 
representing the '^ primitive nebula" of astronomy. Imagine a sphere S drawn in 
the nebula, the radius being a, and the pressure at this surface rr. The matter 
inside S is to form a spherical nebula of finite extent, bounded by a sphere over 
which the pressure is tt, and this matter is to be of a density sufficient to warrant us 
in assuming the gas-equations at every point. The surface S will be continually 
traversed by matter, but this will be of no consequence if the losses and gains 
balance in every respect. The matter outside S must supply the pressure " tt, and 
will also, as was explained in the introduction (§ 3), influence the matter inside S by 

its motion. 

Imagine the matter inside S to be executing a small vibration, and consider two 
extreme hypotheses as to the behaviour of the matter outside S. 

Suppose, in the first place, that the matter outside S is such that it and the 
matter inside S together form a perfect spherical nebula at rest. Then the motion 
of the matter outside S is given by the equations of vibration of such a nebula, and 
the influence of this matter upon that inside S is exactly that required in order to 
enable the matter inside S to execute the vibrations given by the equations of an 
infinite nebula. 
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Suppose, next, that the matter outside S consists mainly of molecules or of masses 
of matter which are describing hyperbolic or parabolic orbits, or which come from 
infinity and after rebounding from the nebula return to infinity. Suppose, further, 
that the interval during which such a mass is appreciably under the influence of the 
matter inside S is so small that it is not appreciably afiected by the motion of the 
latter. In this case the matter outside S may be regarded as arranged at random, 
independently of the vibrations of the matter inside S ; it will not, as under our first 
supposition, take up the motion of the matter inside S to any appreciable extent; 
Hence the matter outside S will exert no force upon that inside S except the 
constant pressure tt, and the vibrations of the matter inside S will be those of a 
spherical nebula of finite size, bounded by a surface at constant pressure tt. 

These two extreme hypotheses lead, as we can now see, to the two conceptions of a 
nebula put forward in § 4. In nature the truth will lie somewhere between these 
two hypotheses, and it is by no means easy to decide which of the two gives the 
better representation of an actual nebula. We shall, however, be within the limits 
of safety if we assert of an actual nebula only those propositions which are true of 
both our ideal nebulae. 

§ 38. We may accordingly sum up as follows :— 

(i.) A nebula at rest and in absolute equilibrium in a spherical configuration will 

always be stable, 
(ii.) Such a nebula may become unstable as soon as the temperature-lag is taken 

into account, 
(iii.) There will be a linear series of configurations of relative equilibrium of a 
rotating nebula, starting from a non-rotating spherical nebula (supposed 
stable), and such tha.t the configuration is symmetrical about the axis of 
rotation. This linear series corresponds to the series of Maclaurin 
spheroids, 
(iv.) The first point of bifurcation on this series occurs for a comparatively small 

value of the angular rotation. 
(v.) The second series through this point is one in which the configurations 
possess only two planes of symmetry. Initially the configuration is such 
that the equations to the surfaces of equal density contain only terms in 
the first harmonic in addition to those required by the angular rotation, 
(vi.) There is a linear series which corresponds to the series of Jacobian ellipsoids, 
each configuration possessing three planes of symmetry. The point of 
bifurcation at which this series meets the series mentioned in (iii.) is a point 
at which the angular rotation is much larger than that at the point of 
bifurcation mentioned in (iv.). 
(vii.) This latter linear series appears to be always unstable. 
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The Unsymmetrical CoNriGuiiATioNS of a Nebula. 

The Second Series of Equilibrium Configurations. 

§ 39. Let us now try to examine the second series of equilibrium configurations, 
which, as we have seen, is a series of stable configurations replacing the series of 
Jacobian ellipsoids. In this way we shall be able to gather some evidence with a 
view to forming a judgment whether the behaviour of the nebula after leaving the 
symmetrical configuration is such as is required by the nebular hypothesis. 

Let us suppose, in the first instance, that the symmetrical configuration from which 
this series starts is one in which there is no rotation, so that the configuration is one 
of perfect spherical symmetry. If the nebula is one in which cooling takes place 
very slowly, the configuration of the nebula will always be very approximately an 
equilibrium configuration. This configuration will be one of the spherically 
symmetrical series until the first point of bifurcation is reached ; after thi^ the 
configuration will change so as to move along the other series, which passes through 
this point. 

Now we have already found the manner in which the configuration first diverges 
from spherical symmetry ; in other words, we have a knowledge of the unsymmetrical 
series in the immediate neighbourhood of the point of bifurcation. If then, we can, 
by some method of continued approximation, obtain a more extended knowledge of 
this series of configurations, we shall be able to trace the motion of a nebula w^hicli 
is cooling with infinite slowness, and in this way form some idea of the motion to be 
expected in the more general case. 

Let us assume, as a general form for the '^ series " now under discussion, 

P ^^ /^o "I Pi^i \ pi^-i'T' Ps^s + • • • • ♦ • ' * (115)5 

where P^ is the zonal harmonic of order ^, and ^q, /)|, p^ '^^^ functions of /' and of 
some parameter a. This parameter determines the position of any particular 
configuration in the series. We shall suppose that at the point of bifurcation a = 0, 
and we then know that Avhen a is very small the limiting form of p is 

In the notation Avhich has been in use throughout the paper, we find that 
corresponding to the density distribution given by equation (72) the gravitational 
potential at the point r^ m 

Yr=e, + e,v, + d,?, + e.j>, + (no), 

where 

VOL. CXCIX.^ — A. a 
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The functions p|, p^ . . . . are to be determined from the condition that V and p 
shall satisfy the three equations of equilibrium, which are of the form 

1 dm clY 

p dx clx 



An Isothermal NehulcL 

I 40. Let us suppose, for the sake of simplicity, that the nebula is at uniform 
temperature, and extends from r = to r = go « We have already seen (equation (77)) 
that the critical vibration for a nebula initially isothermal, is one in which the nebula 
remains isothermal. Hence it follows that if a nebula changes its configuration 
through coming to a point of bifurcation, when moving on a series of isothermal 
and spherical configurations, then the new series will also be one in which the 
equilibrium is isothermal. 

We may now write m = Kp, where k is a constant, and the three equations of 
equilibrium become equivalent to the single equation, 

K log p = V + c , 
or 

Now the series in question is, as we have seen, approximately represented, near to 
the point of bifurcation, by taking only two terms of (115), and consequently only 
two terms of (116). In this case equation (118) becomes : 






K 



l + ^j-p, + il^^^X+-..\ . . . (119). 



Equating coefficients, we find that pQ is given by the equation 



Po - e ^ , 



the same equation as in the case of perfect spherical symmetry. Also p^ is given by 
the equation 



Pi = '"-' e. 



It will be easily verified that this equation is exactly equivalent to our former 
equation (38). The equation contains an arbitrary multipHer in its solution. This 
may be taken to be a, the parameter of the series, so that we may write 



p^ =::: acTj,, 
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where cti is a conapletely determined function of r. Thus, as far as a, the solution is 
seen to be 

P = Po + ctcr^Py 

We shall now show that, as far as a^, the solution is 

/> ^ Po "f ^^^^03 + <^^<^iPi + ^^^p.Ps ...... (120), 

The substitution of this in equation (118) leads to 

= e"?^l, +°*'P, + i(''+!p,y + . . . +'*P, + . . .+?'*» + . ..1, 

where ^^ stands in the same relation to (Ti as does 6^ to pj. The right-hand member 
of this equation is equal to 



l + xi^W + ^. + ^P^ + /^lt.+ x^>^ + ...l 



/<;" fC K ^ \ K ' ^ t^ 



in which the unwritten terms are of degree at least equal to 3 in a. 
Neglecting a^ the equ.ation is satisfied if 



' i^i' , <!> 



03 



Po = e«, o-o3 = po iti-H-f' (121) 



fC IC' ^ ^ 

These equations determine ctq^ and 0*2 uniqu.ely. 

It is obvious that this method is capable of indefinite extension, and that the 
general form of configuration in the series will be given by 

p = Po + ^^ctq^ + ^^<^oi + . . . + {(^(Ti + aVi3 + a^cTjs + . . .) P^ 

+ (aV^ + aVg, + _ . ) P^ + {a^ao^ + a^a-^^^ + . . . ) P3 + &c. . (124). 

§ 41. Let us examine in greater detail the solution as far as a^, this being given by 
equation (115). The important question, as will be seen later, is the determination 
of the sign of or^. We therefore pass at once to the consideration of equation (123). 
Written out in full, this becomes 



"^^"^ bK 1.^ 



Jo Jo ^' 



4-Po [ 1 f '• ^^y, ^^, + r^ f ; ^ ' 
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This equation may be transformed in the same way as equation (39). 

write 

tcr f 



If we 



V 



(T 



Po 



3 



/c 



Q /J 



we find that the above equation is equivalent to (cvf, equations (47)^ (48)^ (49), 
(54), (56)) 

. • • • { L jU *J J , 



dhj 6?/ 
together with the two equations 



Airo-oV 
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o 
OK 



'1 ^ 

r^ dr 



(yr^) 







:Q0 



d^ 

dr 



^O _____ ^y^. 3^ 
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(126), 



(127). 



Writing 



ii 



27rpQr'^ 



K 



equation (119) becomes 



dh/ 

dr'^ 



i (6 - 2^0 



• > 
OK," 



(128). 



Referring to the table of values for ii^ which will be found on p. 15 of 
Professor Darwin's paper, it appears that u increases from a zero value at the origin 
up to a maximum value of about 1"66 ; it then decreases to a minimum of about '8, 
and after this increases to 1, its value at infinity. Thus the factor 6 •— • 2ii has a range 
of values from 6 to about 2f . 

Now the solution of 



d') 



y 



,i'v 



'- n{n+ 1) 



^7>o0r 






(129) 



is easily found to be 

4:7TT 



y 



2n + 1 lr^^+^- Jo 3/c3 



/»00 



dr^ + r''^ 



Po^i' 



J 11 ». 



]fc\^n~l 



dr' 



> + C,r-^ + C^r^^+^ . (130)^ 



in which Cj and C^ are constants of integration, which may at once be put equal to 
zero, if n is positive, and if y is to satisfy conditions (126) and (127). 

Comparing (128) with (129), we see that if u had a constant value t/g at every 
point of the nebula, the value of y would be given by equation (130), in which C^, C^ 
would be put equal to zero, and n would be the positive root of 

n{n -^ 1) =:: Q — 2uq, 



provided only that 6 — 2uq were positive. 
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For the range of values for Hq from iCq = to it^ == 1*66, the vahie of n would have 
a range of values from 2 to 1*2. Thus the form of solution is materially the same for 
all of these values o{ Uq. It will be seen without diiSiculty that the solution of (128), 
in which u has not a constant value, but varies over the range from to 1*66 as r 
varies, will be such that the graph expressing y as a function of r will present the 
same features as are common to the graphs gi^en by equation (130) for ranges of n 
from 2 to 1'2. 

Now the value of y given by equation (130) is positive for all values of r, hence we 
infer that the solution of (128) is such that y is positive for all values of r. We 
therefore have, for all values of r. 



= ^ —J- + ^ positive quantity. 






so that cr^ is positive for all values of r, 

§ 42. We therefore see that the initial motion, in which u and A are each 
proportional to the first harmonic, will first break down owing to the introduction of 
terms involving the second harmonic. The sign of these terms is such that there is, 
in all the shells of which the nebula is composed, a diminution of density in the 
equatorial regions, and a condensation at both poles, which must be added to that 
given by the terms involving the first harmonic. 

The nature of this motion will become clearer upon a reference to fig. 3. This 
figure consists of the four curves^^' 

r = cIq r = aQ + ^iPi 

r = ao + c^iiPi + CI2P9. ^' = ^k + ^ nPi + <^\^2^ 

and these may be supposed to represent curves of equal density in the three stages. 
It is easy to see that of the pear-shaped surfaces of equal density, the equations of 
which contain the two first harmonics, some will be turned in one direction, and some 
in the other. For if they were all turned in the same direction the centre of gravity 
could no longer remain at the centre of co-ordinates. Thus, if the narrow ends of 
these pear-shaped figures point in one direction at infinity, we must, as we go 
inwards, come to a place at which they have the transition shape, namely, ellipsoids 
of revolution, and after this they will point in the opposite direction. 

It appears, therefore, that the initial motion is such as to suggest the ultimate 
division of the nebula into two parts, this division being effected by the outer layers 
condensing about one radius of the nebula, so as to leave room for the ejection of a 

■^" The particular values for which the curves are drawn are in the ratio <:^^^ = 11, ^i = 2, ctn = 5, cio = 2, 
a'li = 7, a' 2 = 4. Thus the equation of the last curves are in polar co-ordinates, 

r r=: 1^ (10 + 5 COS ^ + 3 cos2 0^^^ ^. ^ ^ (9 + 7 cos 6* 4- Q cos2 0). 
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central niiclens in the direction of the opposite radius. Whether or not actual 
separation takes place would probably depend on the amount of the angular velocity. 

It is of interest to compare the result just arrived at, with the corresponding 
result found by Potncaiie for the motion when an ellipsoid of Jacobi first becomes 
unstable.^^ This is described as follows :— 

" La plus grande portion de la matiere semble se rapprocher de la forme spherique, 
tandis que la plus petite portion de cette meme matiere sort de rellipsoide par 
Textremite du grand axe, comme si elle voulait se separer de la masse principale,^^ 

Thus, although the initial motions are, since they start from different configurations, 
necessarily different, yet it would seem as if the final result was very much the same 
in the two cases. In either case we have a diminution of matter in the equatorial 
regions, suggesting the ultimate division of the mass into two, and in each case these 




Fi 



g. 3, 



two masses are of unequal size, a result which could hardly have been foreseen 
without analysis. 

§ 43. If the rate of cooling of a nebula is appreciable, the motion will not be along 
a ^^ series" of equilibrium configurations. The value of ^9, the frequency which is 
nearest to instability, will be changing at a finite rate, and may run to some distance 
beyond the zero value, before the deviation of the nebula from the spherical shape is 
sufficient to invalidate the analysis of our paper« In this case we can imagine the 
first unstable vibration, that for which 2^ = 0, being overtaken by other unstable 
vibrations of greater and greater frequency, the corresponding velocity of divergence 



^ i 



Acta Mathematica/ vol. 7, p. 347. 
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from spherical symmetry becoming coiitimially greater. It is therefore quite 
conceivable that the motion may become adiabatic at an early stage, and it is possible 
that it may be better imagined as a collapse or explosion, rather than as a gradual 
slipping from a spherical state of equilibrium into and through a series of 
unsymmetrical states of equilibrium. 

But an examination of the physical character of the motion will show that in this 
extreme case, as also in any intermediate case, the motion must be, in its essentials, 
the same as that which has been found for the other extreme case, namely, that of 
infinitely slow cooling and perfect thermal equilibrium. In the spherical state, the 
outermost layers of gas may be regarded as stretched out in opposition to their 
gravitational attractions, being maintained in this state by the elasticity of the gas. 
The balance between these two agencies (which is, speaking loosely, measured by the 
stability function, u^) must be supposed to be continually changing, and instability 
always results from the same cause, namely, that the elasticity of these outer layers 
becomes inadequate to resist the gravitational tendency to collapse. In every case 
the outer layei'S concentrate about a single radius of the nebula, the axis of harmonics 
(^ = in equation (72)) and so increase the pressure along this radius, while 
decreasing that along the opposite radius {d = it). This pressure acting upon the 
inner layers of gas and the core sets them in motion, and in this way we have the 
tendency to separation into two nebula. 

A Nebula in ^' Isothermal-acUahatic ^' JEquilihrmm. 

§ 44. A nebula which consists of an isothermal nucleus with a layer in convective 
equilibrium above it, is said to be in '' isothermal-adiabatic '' equihbrium. At the 
surface at which the law changes from the adiabatic to the isothermal, the quantities 
m, T and p must all be continuous. 

The isothermal part is capable of executing a vibration of frequency p = while 
remaining in isothermal equilibrium throughout, provided the forces acting upon it 
from the adiabatic part are* the same as would act if the adiabatic part were replaced 
by an isothermal part in such a way that the whole made up an infinite isothermal 
nebula. If the nebula is rotating, the amplitude of vibration of the infinite nebula 
will vanish at infinity proportionally to some inverse power of n this power increasing 
with the rotation. For sufficiently large rotations, the vibrations may be regarded 
as inappreciable except over the original isothermal nucleus, so that the vibration is 
approximately unaltered when the outer layers are again replaced by layers in 
convective equilibrium. 

We see, therefore, that an ^' isothermal-adiabatic " nebula may become unstable, for 
sufficiently large rotations, through a vibration of order n := I. No attempt is made 
to obtain any numerical results. We can, however, follow up the subsequent motion 
in the same way as in the case of an isothermal nebula. 
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Over the part of the nebula which is in adiabatic equiUbrium, the relation between 
density and pressure is 

where c is a constant, so that the equations of equilibrium become 

__^ dp dV ., 
" dx dx 

and are therefore equivalent to the single equation 
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where Vq is the potential of the outer boundary of the nebula. This takes the form 
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It is obvious that equation (124) again gives the general form of solution, and that, 
as far as a^, the equations are {cf, equations I2I3 123) 
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Writing k for cyp(/ ', vfG see that equations (132) and (133) may be written 
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These are equations similar to (122) and (123); the last term in (136) is difierent 
from the last term in (123), but both terms agree in being invaluably positive. Hence 
it appears that the question of the sign of or^ turns, as in § 37, upon the sign of the 
factor (6 — 2u). We can no longer actually evaluate this factor, as in § 37, but it 
seems to be safe to infer from analogy that it will be positive at every point, and this 
in turn shows that cr^ must be positive at every point. Hence it appears probable 
that the motion will be that described in § 38. 

Rotating Nebula, 

§ 45. The equations of an unsymmetrical series starting from a symmetrical 
configuration in which there is a finite amount of rotation would be extremely 
complicated, and no attempt to handle them is made in this paper. The correction 
for a small rotation will clearly consist merely of an increase in the terms containing 
the second harmonic, so that the general shape of the curves will be similar to that 
of the last two curves of fis^. 3. 

Little difficulty will be experienced in imagining the shape of curves appropriate 
to larger rotations. 

PiiOBLEMs OF Cosmic Evolution. 

Infinite Space filled tvith Hatter, 

§ 46. A limiting solution of the equations of equilibrium (corresponding to A = oo , 
B = CO in equation (114)) gives a nebula in which the density is constant every- 
where. This solution may be supposed to represent infinite space filled with matter 
distributed at random. If space has no boundary there is presumably no need to 
satisfy a boundary-equation at infinity, so that p may have any value ; if, however, 
this equation must be satisfied the only solution is p = 0. 

Let us consider the former case. Space is filled with a medium of mean density p 
and of mean temperature T. Since the space under consideration is infinite, we may 
measure linear distances on any scale we please, and, by taking this scale sufficiently 
great, we can cause all irregularities in density and temperature to disappear. We 
may, therefore, suppose at once that the density and temperature have the constant 
values p and T. 

The equations of motion for small displacements referred to rectangular axes are, 
in the old notation (of, § 6), since V^ and m^ are constants, 



^^ _ ^V^ _ 1 dm^ 
cW dx pQ dx 

or, operating with djdx, d.jdy, djdz, and adding 
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^ -^ VW - W . . . . . . . . {138). 

Since W is the gravitational potential of a distribution of density — Ap {cf, § 6), we 

have 

VW^ -_- 47i'Ap5 ...... o . . (139), 

while if we suppose, for the sake of simplicity, that the motion is adiabatic, so that 
the ratio of pressure to density changes at a constant rate /c, we have {cf, equation 

(3), p. 5) 

V' ^T ^^^^ kS/^p "~" ~"~ '^'^ ^ 

Hence equation (138) becomes 
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The simplest solution of this is of the form 
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and the general solution can be built up by superposition of such solutions. 

Now solution (141) gives A = at infinity, provided q is real, and therefore 
provided p^ + ^Trp is positive, a condition, which admits of p being imaginary. 
There is therefore a possible motion, which consists of a concentration of matter 
about some point, the amount of this concentration vanishing at infinity, and the 
amount at any point increasing, in the initial stages, exponentially with the time. 

We conclude, therefore, that a uniform distribution in space will be unstable, 
independently of the mean temperature or density of this distribution. '^' 

The Evolution of NehulcB. 

§ 47. We can also see that a distribution of matter which is symmetrical about a 
single point willbe equally unstable. For, if this distribution of matter were perfectly 

■* An interesting field of speculation is opened by regarding the stars themselves as molecules of a 
quasi-gas. If space were Euclidean and unbounded, there would be no objection to this procedure, and 
we should be led to the conclusion that the matter of the universe must become more and more concen- 
trated in the course of time. If space is non-Euclidean, this concentration might reach a limit as soon as 
the coarsegrainedness of the structure attained a value so great that the distance between individual 
units became comparable with' the radii of curvature of space. In any case, it may reach a limit as soon 
as an appreciable fraction of the space in question becomes occupied by matter. 
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homogeneous, the whole mass of matter would form a spherical nebula of literally 
infinite extent, and would therefore be in neutral equilibrium. The introduction of 
even the smallest irregularities into this structure is equivalent to the application ot 
an external field of force. This, as has already been seen, will destroy the spherical 
symmetry, and it can easily be seen that the motion from spherical symmetry is such 
as to lead to a concentration of matter about points of maximum density. 

It appears, therefore, that the configuration which will naturally be assumed by an 
infinite mass of matter in the gaseous or meteoritic state consists of a number of 
nebulae (ie., clusters round points of maximum density). We may either suppose 
the outer regions of these nebulae to overlap, each nebula satisfying the gas -equations 
by being of infinite extent, or we may suppose the nebulae to be distinct and of finite 
size, the interstices being filled by meteorites or other matter, which by continual 
bombardment upon the surfaces of the nebulae supply the pressure which is required 
at these surfaces by the equations of equilibrium. 

§ 48. What, we may inquire, will determine the linear scale upon ¥/hich these 
nebulae are formed ? Three quantities only can be concerned : y the gravitational 
constant, p the mean density, and XT the mean (plasticity. Now these quantities 
can combine in only one way so as to form a length, namely, through the expression 

/xT 

of which the dimensions will be readily verified to be unity in length, and zero in 
mass and time. We conclude, then, that the distance between adjacent nebulae will 
be comparable with the above expression. 

Now the value of y is 65 X 10~^, and if we assume the primitive temperature 
to be comparable with 1000° (absolute) we may take XT = 10^ (corresponding 
accurately to an absolute temperature of 350° for air, 2800° for hydrogen). If we 
take the sun s diameter as a temporary unit of length, the earth's orbit is (roughly) 
of diameter 200. If we suppose the fixed stars to be at an average parallactic 
distance of 0*5^^ apart, measured with respect to the earth's orbit, we find for their 
mean distance apart, about 4 X 10*^ sun's radii„ The density of the sun being, 
in C.G.S. units, roughly equal to unity, we may, to the best of our knowledge, 
suppose the mean density of the primitive distribution of matter to be about 
(4 X 10*^)"^, or say 10"®^. Substituting these values for y, XT and p, we find as the 
scale of length a quantity of the order ot 10^^'*^ centims. The distance which 
corresponds to a parallax of 0*5^'' would be about 10^^*^' centims. It will therefore be 
seen that we are dealing with distances which are of the astronomical order of 
magnitude. 
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The Evolution of Planelary Systems, 

§ 49. Let us now regard a single centre, together with the matter collected round 
it, as the spherical nebula which is the subject of discussion. On account of the 
way in which it has been formed, this nebula will, in general, be endowed with 
a certain amount of angular momentum, ¥/e have seen that a primitive nebula of 
this kind may be supposed, under certain conditions^ to become unstable. We have 
also seen that the motion, when the nebula becomes unstable, is such as to strongly 
suggest the ejection of a satellite. 

As a nebula cools the rotation increases, owing to the contraction of the nebula, 
and Xl also increases. Thus the quantity oysX^'P, which measures the rotational 
tendency to instability, has a double cause of increase ; firstly owing to the increase 
in O, and secondly owing to the decrease in T. We can accordingly imagine the 
primitive nebula becoming unstable time after time, throwing off a satellite each time. 

In the usually accepted form of the nebular hypothesis, the rotation is supposed to 
be the sole cause of instability, so that the system resulting from a single nebula 
ought theoretically to be entirely symmetrica] about an axis. On the view of the 
present paper, there is no reason for expecting this symmetry. For large rotations 
of the primitive nebula, the configuration of the resultant planetary system will 
approximate to perfect symmetry, but for small rotations, a slight irregularity 
occurring at the critical moment, at a point out of the equatorial plane, may produce 
a satellite of which the orbit is far removed from the equatorial plane. 

In conclusion, two particular cases of ^^ irregxilarities " may be referred to. If the 
nebula is penetrated by a wandering meteorite, at a moment at which it is close to a 
state of instability, the presence of the meteorite will constitute an irregularity, and 
may easily result in the formation of a satellite^ And if a quasi- tide is raised in the 
nebula by the presence of a distant mass, the same result may be produced. In the 
former case, the plane of the satellite would, if the rotation is sufiiciently small, 
be largely determined by the path of the meteorite ; in the second case, by the 
position (or path) of the attracting mass. It would not, in either case^ depend much 
upon the axis of rotation of the nebula. 

Conclusion. 

§ 50. To sum up, it appears that the behaviour of a gaseous nebula differs in 
at least two important respects from that of an incompressible liquid. In the first 
place, it differs as regards the amount of rotation which is required to produce 
instability, and, in the second place, it differs as regards the disposition of the orbits 
of the planets which will be formed out of the primitive nebula. It will be noticed 
that no definite numerical results have been obtained ; my aim has been to obtain 
qualitative rather than quantitative results, so as to show, if possible, that the 
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results to be expected for a gaseous nebula are of so much more general a kind than 
those usually inferred from the analogy of a liquid mass, that no difficulty need be 
experienced in referring existent planetary systems to a nebular or meteoritic origin, 
on the ground that the configurations of these systems are not such as could have 
originated out of a rotating mass of liquid. 

In conclusion, I wish to express my indebtedness to Professor Darwin for much 
assistance which 1 have received from him throughout the com^se of my work^ 
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